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ABSTRACT 


2-44-61-1 


This  study  is  concerned  with  the  question  of:  (1)  The  feasibility  of  inject¬ 
ing  1-Mev  electrons  into  the  geomagnetic  field  from  a  satellite-borne  electron 
accelerator,  (2)  che  types  of  geomagnetic  experiments  which  can  be  performed 
with  the  injected  electron  beam,  and  (3)  the  import  of  the  experiments. 

In  summary,  the  results  of  the  study  are  as  follows:  (1)  A  satellite  will 
require  a  payload  capability  of  about  1,  000  pounds  in  order  to  carry  an  elec¬ 
tron  accelerator  with  the  desired  characteristics.  (2)  Two  different  types  of 
accelerators  are  required,  each  for  specific  experiments.  There  are  com¬ 
mercial  accelerators  which  would  be  suitable  for  the  injection  experiments. 
One  is  a  1-Mev  pulsed  accelerator,  100  amps  for  4  x  10‘^  sec  with  a  repeti¬ 
tion  rate  of  about  one  per  second.  The  other  is  a  1-Mev  d-c  accelerator  with 

_3 

an  output  current  of  10  amps.  (3)  The  electron  injection  and  trapping  will 
not  be  seriously  affected  by  either  satellite  charging  or  the  diamagnetic  field 
of  the  injected  electrons.  (4)  At  altitudes  greater  than  about  four  earth  radii, 
magnetic  shells  can  be  filled  with  trapped  electrons,  and  their  trapping  life¬ 
time  can  be  subsequently  measured.  (5)  The  bounce  time  of  trapped  electrons 
can  be  measured  to  an  accuracy  of  better  than  1  percent.  (6)  The  geographic 
position  of  magnetic  conjugate  points  can  be  measured,  but  the  accuracy  of 
the  measurements  is  dependent  on  the  type  of  radar  system  employed  in  the 
experiment.  With  either  aircraft-borne  radar  or  a  complex  ground  radar 
system,  it  sho'.ld  be  possible  to  measure  the  position  to  conjugate  point  to 
within  10  km. 
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EXPLANATION  OF  UNITS 


Throughout  this  report,  the  mks  system  of  units  has  been  used.  In  a  few  cases,  for 
the  purpose  of  conventionality,  numerical  values  of  quantities  are  listed  in  cgs  units. 
When  these  quantities  appear  in  equations,  however,  the  mks  values  are  to  be  used. 

A  possible  source  of  confusion  is  the  assignment  of  two  magnetic  moments.  One, 

designated  as  m  ,  is  the  moment  of  a  permanent-magnet  dipole  and  has  units  of 

webers-meter.  The  other,  designated  as  M  ,  is  the  moment  of  a  current  loop  and 

2 

has  units  of  ampere-meter  .  The  two  moments  are  related  by  the  equation 

m  =  Mq  M  (1) 

where  is  the  vacuum  permeability. 
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Section  1 
INTRODUCTION 


The  purpose  of  this  study  is  to  examine  the  feasibility  of  Introducing  a  measurable  number 
of  high  energy  electrons  (the  order  of  1  Mev)  Into  the  region  of  the  earth's  magneUc  field 
above  the  significant  atmosphere,  and  to  determine  the  types  of  experiments  which  can  be 
performed  using  these  electrons  as  a  rejearch  tool.  The  Argus  experiment  has  demon¬ 
strated  that  It  Is  possible  to  Inject  artificially  electrons  Into  the  earth's  field.  This  study, 
however.  Is  concerned  with  InJecUon  mechanisms  otlier  than  nuclear  explosions  and  has 
been  directed  toward  more  precise  experiments  In  which  the  position  of  Injection,  the 
electron  energy,  the  number  of  electrons,  and  the  Injection  pitch-angle  can  be  controlled. 


The  work  Is  divided  Into  four  parts.  First,  for  the  purpose  of  definition,  a  general 
description  will  be  given  of  the  facts  known  about  the  earth's  magnetic  field,  the  solar 
Interaction  with  the  field,  charged  particle  motion  In  the  field,  and  the  Van  Allen  belt. 
Second,  a  compilation  of  techniques  for  Injecting  the  electrons  Into  the  field  Is  given  with 
particular  attention  being  given  to  a  satellite-borne  electron  accelerator.  Third,  an 
Investigation  of  the  experiments  which  can  be  performed  with  the  Injected  electron  beam 
Is  reported  relating  the  results  of  the  experiments  to  the  problems  of  current  Interest  In 
tlie  fields  of  geomagnetism  and  the  magnetic  trapping  of  charged  particles.  Finally,  an 
examination  Is  made  of  those  problems  which  might  limit  the  feasibility  of  electron  Injection 
Involving  the  areas  of  satellite  charging  and  the  effect  of  the  diamagnetic  field  associated 
with  the  injected  current.  In  each  of  the  sections,  the  work  Is  not  Intended  to  be  ex¬ 
haustive,  but  rather  a  thorough  study  of  the  obvious  problems  and  experiments.  As  the 
program  reaches  fruition,  a  more  Intensive  analysis  will  undoubtedly  turn  up  other 
useful  experiments  and  perhaps  unforeseen  problems. 


1-1 


LOCKHEED  AIRCRAFT  CORPORATION 


MISSILES  and  SPACE  DIVISION 


2-4'1-Gl-l 


Section  2 

THE  EARTH'S  MAGNETIC  FIELD* 

The  gross  features  of  the  geomagnetic  field  are  similar  to  those  of  a  uniformly 
magnetized  sphere;  tiUs  dipole  approximation  to  the  field  is  sufficiently  accurate  for 
most  of  the  work  in  this  study.  If  the  earth's  field  is  treated  as  that  from  a  dipole 
situated  at  the  earth's  center,  the  dipole  has  a  moment 

17 

1. 13  X  10  weber-meter 

From  Smythe  (Ref.  1),  a  magnetic  dipole  with  moment  m  has  the  field  components 
given  Iwlow  (see  Fig.  2-1  for  definition  of  coordinates). 


Z 


Fig.  2-1  Definition  of  the  Coordinates  of  a 
Typical  Field  Line  of  a  Magnetic  Dipole 

*The  material  in  this  section  and  the  following  section  is  taken  in  part  from  Satellite 
Environment  Handbook  (Ref.  3) 
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«  _  m  sin 

®r  27  3 

r 

^  ^-jn  cos  >!• 

4jr  3 
r 


(2.1a) 

(2.1b) 


(2.  Ic) 


The  magnitude  of  the  Held  is 

B  =  (bJ  *  +  Blf‘  =  a  ^  (l  *  3  81..^  (2.2) 

A  magnetic  line  of  force  has  the  equation 

r  =  r  cos^  Ip  (2.3a) 

e 

4>  ~  constant  (2.3b) 

when  is  the  distance  from  the  dipole  to  the  point  where  the  field  line  intersects  the 
equatorial  plane  {P  =  0)  .  The  total  strength  of  the  field  along  a  given  line  of  force 
is,  therefore, 


B  = 


m  1  (l  +  3  sln^  Ip) 

47r  3  6  , 

r  cos  Ip 


1/2 


(2.4) 


The  angle  y  between  a  line  of  force  and  the  radius  vector  is  given  by 


tan  y  =  ^  cot  Ip 


(2.5) 


TT 

The  inclination  of  the  field  Is  -  Y 
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the  next  section,  it  will  be  necessary  to  have  the  expression  for  the  radius  of  curvature 
a  field  line 


(2.Gn) 


r 

c 


/  2 
\1  +  3  sin  </</ 

(1  +  sin^  0 ) 


(2.6b) 


A  better  approximation  to  the  geomagnetic  field  is  obtained  if  the  condition  is  dropped 
that  the  dipole  is  situated  at  the  earth's  center.  The  best  agreement  (Ref.  2)  with  the 
real  field  is  obtained  if  the  dipole  is  shifted  436  km  from  the  center  towards  the  point 
15. 6”N.  lat. ,  150. 9*E.  long. ;  the  dipole  is  then  tilted  so  that  its  axis  Intersects 
the  earth's  surface  at  the  two  points,  75. 0” S.  lat. ,  120. 4”  E.  long. ,  and  81. 0”  N.  lat. , 

84. 7*  W.  long.  With  an  eccentric,  tilted  dipole,  it  is  sometimes  necessary  to  differentiate 
between  a  geographic  coordinate  system  and  a  magnetic  coordinate  system.  The 
following  notation  will  be  used: 

Geographic  coordinate  system  (spherical  coordinate  system  with  the  origin  at  the  earth's 
center  and  the  Z-axls collinear with  earth's  spin  axis). 


4'  =  latitude,  measured  from  the  geographic  equator 

0  =  colatitude,  measured  from  the  earth's  spin  axis 

4  =  longitude,  measured  from  Greenwich 

R  =  radial  distance,  measured  from  the  earth's  center 

0 

Rg  =  radius  of  the  earth,  6. 37  x  10  m 

R  =  radial  distance  in  units  of  earth  radii,  measured  from  the  earth's  center 
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magnetic  coordinate  system  (sptierical  eoordinnte  system  witli  tl«e  origin  at  the  eartlVs 
magnetic  dipole  and  the  Z-axis  colincar  wiUi  tlie  dipole  axis). 

ff*  =  latitude,  measurecl  from  the  magnetic  equator 
0  =  eoIaUtude,  measured  from  U»e  magnetic  dipole  axis 

=  longitude,  usually  measured  from  approximately  70*\V.  gcognipidc  longi¬ 
tude;  ean  be  arbitrary,  however 
r  =  radial  distance,  measured  from  the  dipole  center 
r  =  radial  distance  in  units  of  earth  radii 

fhe  conjugate  points  of  a  given  field  line  in  the  simple  ease  of  a  centered  dipole  are 
defined  as  the  two  points  which  satisfy  #  -  .  In  the  displaced  and  tilted  system, 

the  conjugate  points  are  defined  as  the  two  points  of  intersection  of  a  field  line  with  the 
surface  of  the  earth. 

If  the  parameters  of  the  geomagnetic  field,  such  as  the  conjugate  points,  must  be  known 
with  precision.the  dipole  model  is  Inadequate.  A  survey  of  the  field  over  the  earth's  sur 
face  shows  many  irregularities,  .and.  in  pl.ices,  Uiese  local  anomalies  may  constitute  an 
important  fraction  of  the  total  field  (Fig.  2-2).  Several  workers  (Ref.  3)  have  made 
spherical  harmonic  analyses  of  the  earth’s  field  using  a  large  but  finite  number  of  terms 
to  match  the  field  at  the  surface  of  the  earth.  As  is  to  be  expected,  this  analysis  shows 
that  the  effects  of  the  local  anomalies  fall  off  rapidly  with  increasing  altitude;  far  above 
the  earth's  surface,  the  dipole  approximation  is  valid.  Vestine  and  Siblev  (Ref.  4)  h.ave 
calculated  the  conjugate  points  for  a  number  of  lines  using  a  48-term  harmonic  an.alysls. 
Their  results  are  shown  in  Fig.  2-3.  Conjugate  points  are  Indicated  by  pairs  of  dots; 
solid  lines  are  the  loci  of  the  conjugate  points. 


Thus  far,  we  have  considered  a  static  geomagnetic  field.  In  reality,  however,  the  direc¬ 
tion  and  magnitude  of  the  field  are  constantly  fluctuating.  The  transient  variations  in  the 
geomagnetic  field  arise  mainly  from  interaction  between  the  geomagnetic  field  and  the 
ionized  gas  moving  out  from  the  sun.  This  gas  which  flows  radially  outward  is  referred 
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Fig.  2-3  Approximate  Intersections  of  Lines  of  Force  of  the  Geomagnetic  Field 
With  the  Earth's  Surface  in  the  Northern  and  Southern  Hemisphere. 
Conjugate  points  are  indicated  by  dots. 
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to  as  the  solar  wind  (Ref.  5).  Tlio  best  direct  evidence  for  the  existence  of  a  continu¬ 
ous  solar  wind  comes  from  the  comet-tail  observations  of  Biermann  (Ref.  6).  The 
solar  wind  will  push  into  Uie  geomagnetic  field  roughly  to  the  point  where  the  kinetic 
energy  density  of  the  solar  wind  is  equal  to  the  magnetic-field  energy  densify.  By 
equating  these  two  energ>'  densities  and  assuming  that  the  magnetic  field  strengUi  is 
given  by 


B 


(2.7) 


where 

B^  =  the  field  strength  at  the  earth's  surface 
Rg  =  the  radius  of  the  earth 

r  »  the  radial  distance  to  the  point  where  the  field  strength  is  B 


we  obtain 


r 


(2.8) 


where  n  is  the  density  of  solar  wind  particles,  m  is  their  average  mass,  and  v  the 
solar  wind  velocity.  Taking  B  =  0. 35  gauss  (3. 5  x  lO"®  webers/m^),  n  =3x10^ 

3  ®  5 

protons  (and  electrons)  per  m  ,  and  v  =  3  x  10  m/sec  for  quiet  solar  conditions,  we 
find  that  the  solar  wind  will  penetrate  to  r  =  7.7  R„.  During  periods  of  solar  activity, 

^  1  A  O 

appropriate  values  for  the  solar  wind  may  be  as  high  as  n  =  10  m  and 
6  , 

V  -  1.5  X  10  m/sec,  leading  to  a  value  r  =  1.7  R^.  'fhe  solar  wind  blowing 
against  the  geomagnetic  field  will  probably  distort  the  field  into  the  shape  shown  in 
Fig.  2-4  (Ref.  7).  The  front  surface  is  determined  by  the  impact  pressure  of  the  solar 
wind  while  the  long  "tail"  is  closed  by  the  pressure  of  the  transverse  thermal  motion 
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of  the  solar  wind  particles.  In  all  cases,  the  plasma  pressure  is  equal  to  ^ the 

interface.  Clearly,  in  the  region  of  high-altitude  satellite  orbits,  such  a  distorted 
field  will  lead  to  large  diurnal  effects. 

Sudden  increases  in  the  solar  wind  pressure  lead  to  the  phenomenon  of  magnetic  storms. 
The  impact  of  the  strengthened  solar  wind  increases  the  distortion  of  the  geomagnetic 
field,  generates  hydromagnetic  waves  which  propagate  to  the  lower  ionosphere,  and 
injects  charged  particles  into  regions  where  tliey  become  magnetically  trapped.  It  is 
tliese  two  major  magnetic  transients,  the  diurnal  effect  and  the  magnetic  storm,  which 
can  be  studied  with  artificially  injected  electron  currents. 
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Section  3 

CHARGED  PARTICLES  IN  THE  GEOMAGNETIC  FIELD 


3. 1  MOTION  OF  A  CHARGED  PARTICLE  TRAPPED  IN  THE  GEOMAGNETIC  FIELD 

Because  of  interest  in  the  effect  of  the  geomagnetic  field  on  primary  cosmic  rays 
approaching  the  cartli,  Uie  problem  of  charged-particle  motion  in  a  magnetic  di|x)lo 
field  has  been  studied  quite  Uioroughly  by  Stormer  (Ref.  8)  and  Alfven  (Ref.  9).  In 
.solving  for  the  ixirticle  orbits,  a  class  of  solutions  was  found  in  which  the  particles 
arc  magnetically  "trapped."  These  solutions  were  considered  to  be  of  academic 
interest  only  because  no  external  source,  such  as  cosmic  rays,  can  inject  particles 
into  trapped  orbits.  Tlie  discovery  of  the  Van  Allen  belt  made  up  of  trapped  electrons 
and  protons  has  brought  renewed  interest  to  the  analysis  of  trapped,  charged-particle 
dynamics  in  the  earth's  field.  A  complete  description  of  trapped  particle  motion  would 
entail  numerical  solutions  of  the  equations  of  motion.  However,  a  good  approximation 
to  the  exact  motion  can  be  obtained  by  making  use  of  quantities  which  are  nearly  con 
served-the  adiabatic  invariants  of  motion. 


A  particle  trapped  in  a  dipole  field  executes  three  principal  motions; 

(1)  A  particle  moves  in  a  circle  of  radius  a  =  p  /(ZeB)  with  an  angular  fre- 

2  ^  'JL 

quency  w  =  , 

(T  +  mc'^) 


where 

m  =  the  mass  of  the  particle 
T  =  its  kinetic  energy 

=  its  component  of  momentum  perpendicular  to  the  local  magnetic  field 
direction 
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c  =  the  velocity  of  light 

7sC  =  Ojo  particle  charge 

B  =  Oie  magnitude  of  the  local  magnetic  field 

For  Von  Alien  rndlaUon  particles,  n^«  r  .  where  r  Is  the  distance  to  the 
center  of  the  earth  from  the  masnetlc  field  line  about  which  the  particle  Is 
spiraling.  This  circular  motion  Is  properly  enllcd  cyclotron  motion; 

le  Oio  cyclotron  radius  and  Is  the  cyclotron  frequency,  flic  center 
of  the  circle  is  called  the  parUclc's  guiding  center. 


(2)  Tho  guiding  center  of  a  particle  ivlmsc  velocity  vector  makes  an  angle  o 
(the  pitch  angle)  with  tho  magnetic  field  direction  at  a  position  where  the 
field  strength  Is  B  moves  nearly  along  the  field  line  to  a  position  where  the 
field  strength  Is  =  B/sln^  a  .  The  position  at  which  B  =  B  Is 
called  the  mirror  point,  and  a  particle  reaching  this  point  will  be  rrflccted 
back  along  the  field  line.  Thus,  a  particle  which  Is  trapped  In  the  geomagnetic 
field,  besides  ckccuUng  small  scale  cyclotron  moUon,  moves  very  nearly 

along  a  field  line  and  bounces  back  and  forth  between  northern  and  southern 
hemisphere  mirror  points. 

(3)  As  a  trapped  particle  bounces  back  and  forth  between  mirror  points,  its 
guiding  center  slowly  drifts  in  longitude,  positively  charged  particles  drift 
westward,  and  negatively  charged  particles  drift  eastward. 


There  arc  three  invariants  which  describe  the  motion  of  a  charged  particle  trapped  In 

a  static  magnetic  field.  Each  of  the  invariants  will  be  described  briefly  as  they  relate  to 

the  Van  Allen  radiation  bell.  For  a  more  detailed  discussion  with  many  references  to 

previous  work  on  this  subject,  see  Northrop  and  Teller  (Ref.  10).  The  three  Invariants 
arc  as  follows: 


(1)  The  first  invariant  is  the  magnetic  moment  invariant, 
of  a  charged  particle  in  a  magnetic  field  is  given  by 


The  magnetic  moment 
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-  1/2  j^c  pj^/B(T  +  mc“)J  .  As  long  as  Uie  particle  moUon  and  field 

configuration  arc  such  U,at  the  magnetic  field  strcngtl,  changes  only  very 
slightly  over  the  distance  the  particle  moves  In  one  cyclotron  period,  U.e 
ma^etle  moment  tvlll  remain  constant,  l.e. ,  pJ/B  is  a  constant  of  the 
0  on.  (  remains  constant  in  a  static  magnetic  field. ) 

=  invariant 

^  '’ll‘^''"*®™''ll  ‘sthccomponentofthopartlcle 

vo  oelty  parallel  to  the  loeal  magnetic  field,  dl  la  the  element  of  length 
along  the  line  offeree,  and  the  Integral  Is  evaluated  between  mirror  points. 

nv^lant  of  the  motion  places  an  additional  eonstraint  on  the  motion  of 
a  rapM  partlele.  It  a  trapped  parUele  drifts  In  the  geomagnetic  field  so 
that  M  and  J  are  eonstant,  the  partlele  must  drift  In  sueh  a  way  that  It 
eventually  returns  to  the  same  field  line  from  whleh  It  started.  The  con¬ 
stancy  of  M  only  requires  that  be  eonstant;  the  additional  constaney 
o  J  requles  that  the  particle  drift  along  a  well  defined  Integral  Invariant 

oalculat  ons  Is  shown  In  Fig.  3-1.  tt.is  figure,  the  Integral  Invariant 
calculations  of  Jensen  et  al. .  |Ilef.  11)  are  used  to  Indicate  some  paths  or 
races  described  by  the  mirror  points  of  particles  trapped  In  the  geomagnetic 


(3)  The  third  Invariant  is  the  nux  Invariant.  The  flux  Invariant  Is  simply  the 

conservation  of  magnetic  flux  Inside  the  Integral  Invariant  surface  defined  In 
the  paragraph  above. 


*  =  /  B  •  ds 
s 


where  ds  is  an  element  of  area, 
joining  the  magnetic  equator  of  the 


The  integral  is  evaluated  over  any  surface 
earth  to  the  integral  invariant  surface.  The 
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I. 


const;, „c.v  of  rciulrcs  II, at.  If  KcoraoBnolic  fi„l,l  «.oro  slon  lv  lo 

coniract  or  ox,,„n.l  (as  happona  ao,„o  ptaaoa  of  BOo.„„Bnotlc  storms), 

11,0  IntcKrai  lavarlanl  surface  mas.  ol.a„so  Its  s.r.„  .accord, agly  so  .l,a.  .s 
IS  coasorvd.  Durtaj  U..s  .noccss,  M  and  J  also  .aaaata  coaslaat.  Aa 
altc  nato  way  of  sUatlag  the  flu.,  „, variant  Is  to  say  II, o  drift,,,,.  p,,r„clc 
follows  slow  inovemonts  of  the  maKijclIc  field. 

l;x,„-,-„„..,„s  doscrldd  in  lalor  sections  arc  dcstpied  lo  measme  some  of  „.e  pa,-amc,..,.s 
.  ssot  ,,„al  will,  mapped  eleCoa  ,„o,lo„.  mil,  this  In  mind,  we  will  e.aamino  In  a,.,,-,. 
<lcl.,ll  .some  of  the  features  of  the  trapped  particle  motion  ,lesci-lbed  al»ve. 

•'i.I.I  Holince  Time 


Tl»!  Mine  for  a  particle  lo  move  from  one  mirror  point  lo  the 
and  flien  reUirn.  Is  called  the  bounce  time. 


conjugate  mirror  |)olnt, 


(3.1) 


In  this  inf  aral,  HI  Is  measured  along  a  field  line,  L  Is  the  distance  from  the  magnetic 

e,|ii..l.ii  to  the  mirror  point,  and  v, ,  Is  the  component  of  the  particle's  velocity  parallel 
to  Ilif  loc.'d  magnollc  field. 


V(  I  =  V  cos  a 

From  the  conservation  of  the  magnetic  moment, 

V|  ,  =  V  (1  -  E/B 

m ' 


(3.2) 


(3.3) 
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If  ;f  =  0  ,  II, c„  V  ia  oonslant  ami  U,c  boonco  lime  la  a  pare  .llpolo  field  Is  glvea  bv 


^  ni  9  1  /o 

'*'n  =  -r  / 

(1  -  B/B 

nr 


(3. 4) 


B/B  =  (1  •»•  3 


C  . 
cos  tji 


m 


m 


cos  (1+3  sin^  ij,  ) 


2  1,  ,1/2 

m' 


(3.5) 


where  la  Ihe  latllude  of  Ihe  nilrror  imlal  (magaelle  eeor.llaatea  are  used  IhrouBhoul), 

aad  Kldntlfles  the  partleular  line  of  foree  belag  eoasldered.  The  Integrand  la 

(3.4)  Is  singular  at  4.  =  4^  ;  Iwwever,  Ihe  singularity  is  Inlegrablo.  The  Inlcgia- 
tlon  has  been  performed  aumerleally  over  most  of  Ihe  range  of  4, ,  aad  la  the  vielnity 
Of  1}'^  an  approximate  expression  has  been  integrated  analytically. 


ean  be  related  to  the  pitch  angle  at  the  equator  a  bv 

e 

6  , 

9  cos  III 

sln^  a  =  - 

e  -  ...„2  A/2 


(3.6) 


(1+3  sin  111  ) 


(3.7) 


Figure  3-2  sho,vs  S(o^)  as  a  fuaetion  of  sin  p.^  .  The  striking  feature  Is  the  eoaslaaey 
of  S  .  For  rough  ealculatlons.the  bouaee  time  ean  be  approaimated  by 


4  r 


B 


(3.8) 


Putting  r^  into  units  of  earth  radii,  and  taking  v  equal  to  the  veloeity  of  light,  the 
bounce  for  a  high  energy  electron  is 


(6.5  X  10  )  r^  see 


(3.9) 
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“eT:  f-">  »•  > 


3.1.2  Longitudinal  Drift  Time 


ar..  ,0..,,.  ^mchI  penTii'r;:  r:;t 


V  =  ^  ( 1  2  2  \ 

^  ‘^c*’c  y  ^  ''ly 


(3. 10) 


w  Is  the  particle's  cyclotron  angular  frequency  and  r  lo  «i. 


^  =  D 
dt  r  cos  Ip 

The  average  value  of  over  one  bounce  time  Is 


(3.11) 


d£ 

clt 


=  f  /"’«.>  iLllstaiu^ 


B  0 


V|| 


2 

cos  ijl 


(3.12a) 


-  6p 

r_  e  r  jj 
B  e  e 


T„  e  r  B  j 


''m  (1  +  sin^^^)  (2  -  B/B  )  cos^ 

-  m  ^ 


where  p  Is  the  electron's  momentum  and  r  la  fu 

the  eaualn.  ».  ,1,  ,  .  ®e  «>'  magnitude  of  the  magnetic  field 

me  equator  at  the  point  r  =  r  The  drift  time  is 


on 


T  =  -2jl 
D  ^ 

dt 


(3.13a) 
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"D 


4  e  r^B 

4ir  0  e 

3  pv 


/ 


in 

coiLilLlAsisitlL^ 

(1  -  B/B 


/■"dO 

i  0  -  (1  +  3  sin^ 


(3. 13b) 


Welch  and  Whitaker  (Ref.  13)  and  Watao...  et  al.  (Ref.  29)  Indicates  that  the  ratio  of  the 
two  Integrals  In  Eq.  (3. 13b)  is  unity  to  about  .30  percent  over  the  range  of  ,!>  from  0  to  60* 
For  high  energy  electrons,  the  drill  Ume  may.  therefore,  be  approximated  by 


4,  e  r^  B 
T  s  IS.  e  e 

D  3  pc 


(3.14) 


On  the  equator,  the  magnetic  field  Is 


Be  =  (9  X  10^®)  r^^  webers/m^ 


(3. 15) 


For  a  one-Mev  electron,  the  drift  time  is  then 


_  4.2  X  10^ 

- = -  sec 

r 

e 


(3. 16) 


where  is  in  units  of  earth  radii. 
1|000  sec  to  100  sec. 


The  drift  times  have  order  of  magnitude  values  from 


3.1.3  Lifetime 


The  description  of  trapped,  charged-particle  motion,  thus  far,  has  not  Included  any  loss 
mechanisms,  and  the  particles  would  be  trapped  Indefinitely.  It  radiative  energy  loss 
is  negligible,  which  Is  the  case  except  lor  very  high  energy  eleotrons,  the  dominant  loss 
mechanisms  are  scattering  and  Ionisation  In  the  upper  atmosphere.  Since  the  atmospheric 
density  falls  off  nearly  exponentially  with  Increasing  altitude,  most  of  the  energy  loss 
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occurs  when  the  particle  is  nearest  the  earlh,  namely,  at  the  mirror  polnU.  A 
detaUed  doserlpUon  of  the  analysis  of  this  energy  degradation  and  eventual  loss  of  the 
particle  from  a  stable.frapped  orbit  is  beyond  the  scope  of  the  present  study.  The 
enalysis  has  been  made,  however,  by  Walt  and  MacDonald  (Ref.  M).  The  general 
features  of  trapped  parUcle  loss  are: 

(1)  The  parUclc  loses  energy  by  cxeltaUon  of  the  atmospheric  atoms  In 
the  vlolnlbr  of  the  mtrror  points.  The  rate  of  energy  loss  Increases  as 
the  particle's  energy  decreases. 

(2)  The  posiUon  of  the  mirror  points  gradually  decreases  In  alutude,  and 
hence,  goes  into  a  region  of  greater  atmospheric  densiW  and  consequenUy 
greater  energy  loss.  This  lowering  of  the  mirror  points  Is  caused  by 
atmospheric  Coulomb  mulUplo-eeatterlng  over  the  particle's  entire  orbit. 

(3)  There  Is  cataclysmic  parUcle  loss  In  the  vicinib'  of  the 

anomalies,  such  as  the  Cepe  Town  anomaly.  Here,  the  geomagnetic  field 

Is  so  much  weaker  than  at  the  conjugate  mirror  point  that  tho<K  particles 

normally  mirroring  at  relaUvely  low  alUtudes  are  dumped  Into  the  atmosphere 
and  lost. 

(4)  To  demonstrate  the  conservation  of  the  particle's  magnetic  moment,  it 
had  to  be  assumed  that  the  geomagnetic  field  changed  negligibly  over  an 
area  defined  by  one  cyclotron  revolution  of  the  particle.  This  assumption 
breaks  down  at  high  altitudes  where  the  field  is  weak  and  the  qrclotron 
radius  is  large.  If  the  magnetic  moment  is  not  constant,  the  mirror  point 
position  is  not  stable  and  the  particles  would  be  only  quasi-trapped  and 
eventually  would  be  lost  into  the  atmosphere.  It  has  been  suggested  that 
the  outer  altitude  limit  of  the  Van  Allen  electron  belt  is  defined  by  this 
magnetic  moment  non-conservation. 

(5)  During  magnetic  ulorms,  the  earth's  field  Is  compressed,  hydromagnetle 
waves  propagate  along  field  lines,  diamagnetic  ring  currents  ore  generated 
and  the  upper  atmosphere  Is  heated  and  expands,  thereby  Increasing  the 
atmospheric  density  In  the  region  of  the  mirror  points.  Each  of  these  effects 
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can  temporarily  increase  the  loss  of  trapped  particles.  In  particular, 
recent  wrk  (Ref.  15)  has  showm  tliat  the  Itydromagnetlc  waves  can  cause 
the  breakdown  of  the  magnetic  moment  invariance  for  Van  Allen  protons, 
and,  thereby,  set  the  outer  altitude  limit  on  the  proton  bolt. 

3. 2  THE  VAN  ALLEN  BELT 


The  Van  Allen  radiation  bolt  has  been  shown  to  consist  of  two  components:  (1)  a  hard 
proton  component  (the  proton  belt)  centered  at  about  10  x  10^  km  from  the  earth's 
magnetic  axis;  and  (2)  an  electron  component  (tlic  electron  belt)  centered  at  about 
22  >;  10^  km  from  the  earth's  magnetic  axis. 

3. 2. 1  Protons 

Contours  of  constant  proton  flux  are  shown  in  Fig.  3-3.  These  contours  arc  symmet¬ 
rical  about  the  integral  Invariant  equator.  Since  the  geomagnetic  field  is  not  centered 
about  the  earth's  spin  axis,  the  height  of  the  contours  above  the  earth's  surface  is  a 
function  of  longitude.  The  location  of  the  integral  invariant  equator  relative  to  geographic 
coordinates  is  shown  in  Figs.  3-4  and  3-5,  Altitudes  corresponding  to  constant  proton 
flux  at  various  positions  around  the  Integral  Invariant  equator  are  expressed  in  kilom¬ 
eters  by  the  numbers  at  the  arrows.  The  location  of  the  integral  invariant  equator 
and  the  altitudes  were  obtained  from  the  high-order,  spherical -harmonic  analysis 
of  the  geomagnetic  field  carried  out  by  Jensen  et  al.  (Ref.  11).  Figure  3-4  indicates 
the  altitude  at  the  bottom  of  the  radiation  belt  where  the  counting  rate  is  just  above 
cosmic  ray  background.  The  field  strength  here  is  B  =  0.22  gauss.  Similarly, 

Fig.  3-5  indicates  the  altitudes  along  the  integral  invariant  equator  where  the  flux  of 
protons  with  energies  greater  than  40  Mev  is  100  protons  cm"^  sec“^  (about  50  times 
cosmic  ray  background) . 

The  proton  spectrum  has  been  measured  by  Freden  and  White  (Ref.  16)  near  the 
bottom  of  the  belt.  Their  results  are  shown  in  Fig.  3-6.  It  is  expected  that  the  proton 


3-11 


LOCKHEED  AIRCRAFT  CORPORATION 


MISSILES  and  SPACE  DIVISION 


3-12 


LOCKHEED  AIRCRAFT  CORPORATION 


MISSILES  and  SPACE  DIVISION 


L^imioc  fucai 


2-44-61-1 


Fig.  3-4  The  Integral  Invariant  Equator  (B  >  0.22  Gauss).  In 
this  figure  is  shown  the  path  which  a  particle  trapped 
in  the  geomagnetic  field  with  90*  pitch  angle  will  fol¬ 
low.  The  latitude  of  this  path,  called  the  integral  in¬ 
variant  equator,  is  indicated  as  a  function  of  longitude. 
Also  shown  are  the  set  of  altitudes  (km)  where  the 
geomagnetic  field  strength  is  B  >  0.22  Gauss.  This 
field  strength  corresponds  to  the  bottom  edge  of  the 
proton  belt,  where  the  proton  flux  is  barely  detectable 
above  cosmic  ray  background.  The  altitude  shown 
should  be  correct  to  within  ^25  km. 
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east  longitude  (DEG) 

fig.  3-5  The  Integral  Invariant  Equator  /R  -  n  9n  o 

to  Fig.  3-4.  but  the  altltLes  shown  are  fo^a^ff^’  i 

C^uss  At  the  altitudes  indicated  B=0.20 

flux  above  40  Mev  is  lo2  protoS  cm-2  sec-l  Invariant,  the  proton 
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spectrum  will  have  an  aJUtudo  dependence,  and.  consequently,  the  Freden  and 
White  results  will  not  necessarily  hold  throughout  Urn  belt. 

3.2.2  Electrons 


Tho  cleclroi.  belt  conloure  ere  else  shown  in  Fig.  3-3.  These  contours  represent 
on  esUmato  of  the  qulet-div  Hux  ot  electrons  with  energies  greater  than  20  kev. 

There  is,  however,  a  mtgor  uncertain^  In  the  electron  nujt.  The  Dux  values  given 
In  this  figure  are  generally  lower  than  those  given  by  other  authors.  However.  U  more 
conservaUve  values  for  the  electron  holt  are  desired,  the  given  flux  values  should  bo 

niulUpIled  by  10  for  the  outer  zone  and  lo'*  for  the  Inner  rone  as  a  reasonable  upper 
limit  for  the  qulet-daj'  electron  flux. 

The  electron  spectrum  has  been  measured  by  Walt  ct  al.  (Hef.  17)  at  the  "foot"  of 
the  outer  belt.  Tbelr  results  are  shown  In  Fig.  3-7.  As  In  the  proton  case,  these 
results  should  not  be  expected  to  be  valid  throughout  the  belt. 

ObservaUons  have  shown  that  the  proton  belt  Is  stable,  but  the  electron  belt  undergoes 
changes  which  are  associated  with  geomagnetic  storms.  After  the  active  part  of  a 
storm  Is  over  (recover-  phase  -  low  K-Index).  the  counUng  rate  of  Instrumentation 
own  through  the  outer  part  of  the  electron  belt  Increases;  the  magnitude  and  time  of 
occurrence  of  these  variaUons  depend  on  the  ^e  of  detector  used  and  the  characteristics 
0  the  particular  geomagnetic  storm.  Following  the  active  phase  of  the  storm  of 
16  Aug  1959.  the  counters  in  Explorer  VI  showed  a  typical  Increase  of  about  a  factor 
of  ffoe  near  the  middle  of  the  electron  belt.  Another  example  of  possible  Ume  variations 
in  electron  belt  Intensity  Is  shown  by  essenUally  Identical  Instruments  (Anton. 
lype-302  Geiger  counters  with  about  1  gm/cm^  of  shielding)  In  the  space  probe 
Pioneer  fV  and  the  satelUte  Explorer  VI.  The  difference  between  the  counUng  rate  at 

the  peak  of  the  electron  belt  for  Pioneer  IV  (extrapolated)  and  the  pre-storm,  quiet-day 
observations  for  Explorer  VI  is  about  10^ . 
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Fig.  3-7  The  Differential  Energy  l^ectrum  of 
Electrons  at  the  Foot  of  the  Outer  Belt 
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The  Explorer  VI  results  show  U.at  the  Umc  variations  a.  c  very  complex  and  do  not 
fit  Into  any  simple  pattern.  For  example,  during  the  eai-ly  part  of  U.e  storm  of 
16  Aug  1959.  Rosen  et  al.  (Ref.  18)  reported  mi  Increase  In  Uie  outer  r.one  In  the 
counting  rate  of  their  plastic  scintillator  (which  detecta  pai-tlcles  directly  with  high 
efficiency  and  bremsstralilung  wlUi  low  efficiency)  at  U.c  smne  time  that  Fan  et  al. 

(Ref.  19)  and  Ai-noldy  et  al.  (Ref.  2(^  report  a  decrease  In  two  Geiger  counters 
and  one  Ionization  chamber  (which  detect  Van  Allen  electrons  from  Uicir  brcmsstrahlung) . 
Neither  the  cause  for  Uie  relatively  rapid  chamjeM  in  the  observed  intensity  of  the  elec¬ 
tron  belt  nor  the  origin  of  these  electrons  has  been  adequately  e.xplalned. 

The  source  of  the  Van  Allen  belt  will  be  discussed  briefly,  although  it  docs  not  bear 
directly  on  the  matter  of  artificial  injection.  Equations  of  motion  show  that  an 
Individual  charged  particle  cannot  approach  the  dipole  field  from  outside  and  enter  a 
trapped  orbit;  further,  there  are  no  definitive  experiments  pointing  clearly  to  the  source 
of  the  trapped  particles.  There  Is,  however,  good  indirect  evidence  on  the  source  of 
Van  Allen  protons,  namely,  the  beta  decay  of  the  neutron  albedo  caused  by  primary 
cosmic  ray  interacUons  at  the  top  of  the  atmosphere.  Thla  may  also  be  the  source  of 
the  Van  Allen  electrons;  another  likely  source  is  the  solar  wind  plasma.  Although 
individual  particles  cannot  enter  trapped  orbits  from  without,  the  solar  plasma  may 
disrupt  the  geomagnetic  field  locally  and  allow  injection.  The  energy  of  solar  plasma 
electrons  Is  thought  to  be  smaU  compared  with  the  Van  Allen  electrons.  Solar  plasma 
injection,  therefore,  would  have  to  be  accompanied  by  the  acceleration  of  the  newly 
trapped  electrons.  This  could  conceivably  bo  done  by  the  betatron  effect  during  mag¬ 
netic  storms.  The  source  of  the  Van  Allen  electrons  is  one  of  the  major  unanswered 
questions  in  space  physics.  Experiments  with  artificially  injected  electrons,  proposed 
in  the  following  sections,  would  not  result  In  any  direct  evidence  on  this  source.  How¬ 
ever,  an  attempt  would  be  made  to  measure  the  lifetime  of  trapped  electrons.  The 
result  of  the  lifetime  measurement  would  give  an  indication  of  the  required  strength  of 

the  source,  since  a  rough  equilibrium  must  exist  between  the  number  of  electrons 
injected  and  the  number  lost. 
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Section  4 

TECHNIQUES  FOR  ELECTRON  RJJECTION 


,-70 


4. 1  INJECTION  FROM  THE  EARTH'S  SURFACE 

r  “  “  I*"”  »f  "«“»«>■»  from  a  reactor, 

potaw  ky«ard.  would  have  pm-Uclc  energy  and  nu,  Intonal^  great  enough  to  allow  a 

detoclr..lo  number  of  parUolea  to  emerge  from  the  top  of  toe  atmoaphere.  (The  neu¬ 
ron.  ^uld  then  deogy  m,d  be  a  aouree  of  eleetrona.)  An  exrnnioaUon  .  f  toe  number. 

I  volved.  however,  will  ,p.leldy  ^.courage  tola  Idea.  For  eleetrona  In  air,  the  rate 
of  energy  loan  la  2  Mev/gm-om-^  m,d  the  radlaUon  length  la  37  gm/em^*.  If  there 
™re  no  ahower  effeeto.  toe  eleetrona  wmild  have  to  have  an  InlUal  energy  of  2  Dev 
n  order  to  penetrate  toe  atmoaphere  (toe  atmoaphere  la  about  1000  gm/cm^  thick). 

o  ahower  mechaniam  weuld  Increaae  tola  energy  by  roughly  a  factor  of  10.  For 
neutron,  in  air,  toe  mem.  free  path  for  nuclear  InteracUon  ia  about  6  gma/om^  In 
panning  through  toe  atmoaphere,  the  Initial  flux  would  be  degraded  by  a  factor  of  10 

4. 2  EXPLOSIVE  INJECTION  FROM  A  BALLISTIC  ROCKET 

The  A^gua  experiment  him  ahown  that  a  nuclear  exploalon  above  the  atmoaphere  wdl 
depoal  a  meaaurable  number  of  electron.  Into  the  geomagnetic  field.  However,  the 
preaent  bomb  morrtorlum  rule,  out  tola  method.  A  almilar  approach  haa  been 
euggeat^  In  which  a  ahaped-charged  high  exploalve  bomb  containing  an  Intenae  radio- 

and  p^tlally  Ionize  the  3-aource.  The  ionized  component  would  be  trapped  by  the 

earto  a  field;  the  aouree  of  eleetrona  would  be  localized  and  magnetically  aupported 
no  torn  It  would  not  fall  back  Into  the  atmoaphere.  A  rough  calculation  haa  Indicated 
that  the  fraction  of  a  aouree  which  would  be  ionized  by  a  conventional  exploalon  la 
negUgtole  Another  objecUon  to  tola  acheme  (aa  well  aa  to  the  Argua-lype  experiment) 
la  toe  lack  of  control  over  the  number  of  injected  electron,  and  their  energy 
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4. 3  INJECTION  FROM  A  SATELLITE 

Two  methods  of  electron  ii^ecUon  from  a  satellite  have  been  considered.  In  the  first 
method,  ■  radJoacUvo  fl-sourco  would  supply  cicctronsj  no  accelci  uUun  would  bo 
ro<iulred.  The  ojoolod  eloolron  beam  would  be  doflaod  In  direction  and  energy  by 
colllmaUon  and  magnetic  aolcctlon.  This  ,,-stom  la  relaUvcly  simple,  but  rotpilres 
a  source  strength  of  about  100  megacurlcs  for  an  electron-beam  current  of  one  ma. 
(TWs  number  ts  calculated  assuming  that  1  parent  of  the  decay  electrons  wuld  have 
the  appropriate  dlrecUon  and  energy  to  pass  through  the  beam  selector.)  nie  dllfl- 
culttes  Involved  In  handling  such  a  scurce  before  launch,  the  danger  of  contaminaUon 
from  an  abortive  launch,  and  the  problem  of  fallout  when  the  satelitte  evenhially  re¬ 
enters  the  atmosphere  are  so  great  that  this  method  must  be  considered  unfeasible. 

The  second  method  of  satellite  Injection  (and  the  only  injection  scheme  In  which  the 
advantages  outweigh  the  dlsadvantageib  Is  to  use  an  eleetron  aceelerator.  For 
electron  energies  about  1  Mev.  there  are  a  number  of  accelerators  which  are  either 
commercially  available  or  In  the  design  stale.  In  the  experiments  we  propose  In  the 
following  sections,  two  different  types  of  accelerator  will  be  required.  One  should 
be  a  pulsed  type  with  a  high  current  and  low  duty-cycle;  the  other  should  be  a  d-c 
machine.  (The  pulsed  accelerator  cannot  be  used  in  all  situations  because  of  the 
sateUlte-charglng  problem.)  Characteristics  of  two  acceptable  accelerators  arc 
listed  In  Table  4-1.  other  comparable  accelerators  can  undoubtedly  be  found. 

In  the  pulsed  accelerator,  the  Instantaneous  power  dissipation  during  the  pulse  Is 

100  Mw.  This  energy  Is  capacitor-stored,  and  the  weight  of  the  capacitors  and  their 
charging  circuit  is  included  in  the  listed  weight  of  250  lbs. 

The  electrostatic  accelerator  is  made  up  of  a  number  of  Wimshurst-type  electrostatic 
generators  connected  in  series,  rather  than  the  more  conventional  Van  de  Graaff  sys¬ 
tem.  In  designing  a  satellite  to  carry  such  an  accelerator,  allowance  must  be  made 
for  the  gyroscopic  effect  of  the  moving  parts. 
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CHARACTERISTICS  OF  TWO  ELECTRON  ACCELERATORS 


CharactorlsUcs 

Varlan  Pulse- Transformer 
Electron  Accelerator 

High  Voltage  Engineering 
Electrostatic  Accelerator 

Electron  Energy 

1  Mov 

1-2  Mev 

E 

3-4  % 

1% 

Beam  Current 

100  amps  (vlurlng  pulse) 

10  ^  amps 

Rilse  Duration 

4  X  10  ^  sec 

(i.c. 

Repetition  Rate 

0— 2.5/scc 

(1.  c. 

Beam  Size 

,  2 

1  cm 

1  2 

1  cm 

Power  Required 

1. 5  kw 

1. 5  kw 

Weight  (without  power 

supply) 

250  lbs 

COO  lbs 

Size 

3  ft  (Ham.  by  5  ft  length 

approximately  the  same 
as  the  Varlan  accelerator 

Avallablllfy 

1  year 

1  year 

The  power  tor  these  accelerators  will  presumably  be  supplied  by  batteries  Intermlt- 
teutly  charged  by  banka  of  solar  cells.  To  estimate  the  weight  of  the  batteries  and 
the  size  of  the  solar  cell  array  required.  It  will  be  assumed  that  the  aecelerator 
operates  continually,  that  the  batteries  eaa  supply  toe  required  power  for  24  hours 
without  recharging,  and  that  the  solar  cells  will  deliver  n  new  24-hour  charge  In 
a  period  of  12  hours.  The  pertinent  numbers  are  approximately; 

(1)  Rechargeable  batteries  deliver  1  kw-hr  per  15  lbs  of  battery  weight. 

(2)  In  the  spectrum  Interval  over  which  solar  cells  are  sensitive,  the  sun 
supplies  1400  watts/m  above  the  atmosphere. 

(3)  The  solar  cell  energy-conversion  efficiency  is  10  percent. 

(4)  A  solar  cell  bank  weighs  1-1/2  Ibs/ft^. 
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weights  can  probably  be  reduced  by  6o  percent. 

To  oMy  U.0  feaslbllfty  of  varloua  oxporlmcaea  uoIok  U,o  cloctron  bean,  ojoctol  from  a 

afofUfo  w.  wn,  b...  .0  aaa,g„  .peo.fto  cbaraCerlsHoa  b,  «,o  aafoll.u,  and  U.o  o  J- 
tron  accelerator.  It  will  be  nssumed  that 

(1)  For  pulsed  eleetroa  experiments.  U,e  Varlan  leeeloralor  will  be  used: 

100  amps  for  4  pace  with  an  electron  energy  of  1  Mov. 

For  d  e  injeeuon.  the  High  Voltage  Engineering  electrostatle  aeoelemtor 
Will  be  used;  1  nia  at  1  Mcv. 

(2)  TOe  satellite  ean  be  put  Into  any  desiretl  orbit,  and  Its  poslUon  as  a  function 
of  time  is  known, 

(.3)  The  satellite  Is  lnstru.nentcd  to  measure  and  reeonl  Ihe  local  direction  and 
magnitude  of  the  carlh's  magnetic  field. 

(1)  The  aoeelerator  can  l«  oriented  by  ground  control  or  by  a  programmer  so 
that  the  ejected  electron  beam  will  have  a  known  but  variable  angle  with 
respect  to  the  local  direction  of  the  earUi's  field. 

(5)  In  transmuting  experimental  data  from  the  satellite  to  the  ground,  there  Is 
no  limitation  set  by  telemetry  bandwidth  on  the  time  duration  of  signals. 

Expjlments  which  make  use  of  pulsed  Injection  could  ooneelvably  be  carried  out  by 

a  ballistic  .ooket-borneaeoe, orator,  rather  than  from  a  satol.lte,  There  would  be, 

‘■eq^irements  of  a  recovery  capability  and  a  great  number  of 

Je  11^'  ^vantage  to  a  ballistic  rocket  system, 

we  will  restrict  ourselves  to  satellite  injection. 
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Section  5 

DESCRIPTION  OF  THE  EXPERIMENTS 
TO  BE  PERFORMED  WITH  ARTIFICIALLY  INJECTED  ELECTRONS 

5.1  INTRODUCTION 

In  the  following  three  sections,  experiments  to  be  performed  with  an  injected  electron 
beam  are  studied  in  detail.  First,  however,  a  general  descripUon  will  be  given  of  the 
experiments  and  the  types  of  results  which  might  be  expected. 

5.2  FILLED  SHELL  EXPERIMENTS 

An  orbit  is  chosen  such  that  during  part  of  one  revolution  a  satellite  would  be  nearly 
tangent  to  an  integral  invariant  surface  (Section  3. 1).  Electrons  are  injected  during 
this  time  of  tangency  with  an  appropriate  pitch  angle  so  that  they  become  trapped. 
Because  of  the  longitudinal  drift  of  the  electrons,  a  thin  electron  shell  will  spread 
around  the  earth  and  have  a  shape  dictated  by  the  Integral  invariant  surface.  After  the 
shell  has  been  filled,  the  satellite  will  stop  injecting  and  bn  successive  revolutions  will 
pass  repeatedly  through  the  shell.  Using  radiation  detection  instruments  on  board  the 
satellite ,  one  could  then  measure: 

e  The  longitudinal  drift  time,  which  would  be  Inferred  from  the  time 
elapsed  between  the  Initiation  of  Injection  and  first  detection 
e  The  change  of  the  electron-energy  spectrum  with  time  (lifetime) 

•  The  diffusion  of  the  well-defined  edges  of  the  shell 

•  A  diurnal  compression  and  expansion  of  the  field  which  would  shift  the 
position  of  the  shell 

•  The  effect  of  a  magnetic  storm  on  the  shape  of  an  integral  invariant 
surface 

•  The  effect  of  a  magnetic  storm  on  the  trapped  electron-energy  spectrum 
and  the  possible  catastrophic  dumping  of  the  trapped  electrons 
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5.3  bounce-time  MEASUREMENT 


Frra  a  “Millie,  a  pulaa  of  eloolrone  would  be  injected  perpendicular  to  the 

^recUM  of  the  local  nutgnetlc  field.  Theae  electron,  would  mce  along  a  neld  line  to 
tto  con  Jugate  mirror  point  In  the  other  hemlephere.  mul  then  relun.  nearly  to  d..  po.l- 

tte  eatelllte.  the  electron,  (after  one  bounce)  would  mirror  at  a  point  ...oral  kilo¬ 
meter.  fr„„,  u,e  ..telllle.  TO.  would  rule  out  direct  delecUon.  The  proamtee  of  Uto 
electron,  at  the  mirror  point,  however,  could  be  detected  from  their  eynchrolron  nuUa- 

In”’!  «•« 

lime  elapaed  between  Injection  and  detecUon.  A  mcauremenl  of  the  bounce  time  will 

gIVGS 


e  The  average  ,.|„e  of  the  geomagnetic  field  along  a  particular  field  line 
^ch  could  be  compared  to  the  predlcHon.  of  eeveral  geomagneUc  model, 
e  The  maximum  latitude  at  which  field  line.  cIom;  injeoUon  at  higher  latl- 
tudes  would  lead  to  Infinite  bounce  times 

.  The  diurnal  dlatortlon  of  the  field  by  the  imlar  plaema  manlfeatlng  llaelf 
as  a  diurnal  change  In  bounce  time  at  a  constant  latitude 
•  The  effect  on  all  of  the  above  geomagnetic  features  due  to  a  magneUc 


5.4  CONJUGATE-POINT  MEASUREMENT 


r“““'  *  ‘o 

ireollon  of  the  local  magnetic  field.  Theae  elecl^n.  would  move  along  a  field  line  to 
the  conjugate  of  the  Injection  point  In  the  other  hemlapheie.  but  would  not  mirror  be- 

ZeM  I?"  ‘d  r““" 

0  lo  lea, ion  would  be  located  by  radar  renectlon  from  a  ground  atatlon.  Prom  the 
poeltlon  of  the  eatelllte  at  injection  and  the  location  of  the  point  where  the  electron, 
lonlae  the  atmoaphere,  the  conjugate  point,  of  a  particular  field  line  can  be  Inferred 
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As  in  the  provlonsnxporlinon.,,  repented  mensurements  nf  the  position  of  the  eonJoRnte 

points  eould  then  be  eorrolated  with  diumni  nnd  .nngnetle  storm  effeets.  A  proelso 

locnuon  of  eonjngnto  points  Is  not  only  of  selentlfle  Interest,  but  hns  many  potential 

military  appileatlons  sufnelently  Important  to  Justify  the  Implementation  of  the  expert- 
mental  program. 
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Section  G 

% 

THE  FILLING  OF  MAGNETIC  SHELLS 


6.1  INTRODUCTION 

In  this  section,  wo  will  examine  the  feasibility  of  Injecting  and  trapping  a  detectable 
number  of  electrons  In  a  thin  shell  centered  about  an  Integral  Invariant  surface.  It 
will  be  assumed  tliat  the  Injection  pitch-angle  can  bo  controlled  from  the  satellite  so 
that  the  electrons  will  always  be  trapped.  Thus, the  principal  iwlnts  to  be  considered 
are  the  flux  Intensity  which  can  be  achieved  In  a  wcll-dcfinod,  Isolated, thin  shell  and 
the  subsequent  history  of  Intersection  of  the  shell  by  the  satellite.  Throughout  the 
calculation,  the  dipole  model  of  the  geomagnetic  field  will  be  used. 

The  shape  of  the  shell  to  be  filled  will  be  defined  arbitrarily  as  the  volume  lying  be¬ 
tween  the  two  magnetic  Invariant  surfaces  shown  In  Fig.  6-1.  The  outer  surface 


MAGNETIC 

POLE 


crosses  the  magnetic  equator  at  a  distance  from  the  center  of  the  earth  of 
(r^  1-25  a^),  and  the  inner  surface  at  (r^  -  1.25  a^).  The  quantity  a^  is  the 
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cyclotron  radius  of  a  1-Mev  electron  at  the  equatorial  point  (r  =  r  .  =  0)  with  a 

90*  pitch  angle.  ® 


a  =-V- 

c  e  B 

e 

In  order  that  the  satellite  inject  electrons  only  Into  this  shell,  the  Injection  must  be 
limited  to  the  time  the  satellite  spends  in  the  volume  defined  by  the  two  field  lines 
which  cross  the  magnetic  equator  at  (r^  +  0.25  a^)  and  (r  -  0.25  a  ).  The  time 
the  satellite  spends  in  this  volume  is  called  injection  time  T,  and  it  is  the  purpose 
of  the  calculations  which  follow  to  detei  mine  Tj  for  various  types  of  satellite  orbits. 
The  injection  time  is  then  used  to  calculate  the  particle  flux  which  can  be  established 
with  a  given  Injection  current,  and  the  longitudinal  extent  of  the  shell  is  determlncti 
by  comparison  of  the  injection  time  Tj  with  the  electron  longitudinal  drift  time  T  . 

A  filled  shell  is  defined  as  one  with  no  longitudinal  gaps,  so  that  electrons  will  be*^ 
encountered  on  any  subsequent  satellite  passage  through  the  shell. 

The  coordinate  eyetem  shown  In  Fig.  6-2  I.  used  for  the  ealeulatlon.  The  orthogonal 
axes  labelled  with  capital  letters  are  fixed  and  centered  with  respect  to  the  earth, 
with  the  Z  axis  coincident  with  the  earth's  spin  axis  and  the  X  axis  passing  through 
the  Greenwich  meridian.  The  second  set  of  axes  labelled  with  lower  case  letters  Is 
also  fixed  with  respect  to  the  earth,  but  with  origin  at  the  displaced  magnetic  dipole 
position  and  with  the  z  axis  along  the  dipole  moment.  The  satellite  orbits  are  thus 
fixed  with  respect  to  the  earth  system  except  for  rotation  about  the  Z  axis  once  per 
day  at  an  angular  velocity  i  The  dipole  field  Is  fixed  In  both  systems. 

It  is  desired  to  describe  the  satellite  position  P  In  terms  of  the  dipole  radius 

vector,  t .  in  order  to  establish  the  satellite  axlmuthal  range  A»  during  which 

le  satellite  Is  In  the  desired  Injection  region.  The  Injection  time  T,  Is  the  time  to 
cover  this  rnnero  I 


Tr  = 


sat. 


$ 

earth 


where  $ 

sat. 


is  the  angular  velocity  of  the  satellite. 


(6.1) 
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The  expression  for  r  in  the  geographic  coordinate  system  is 

I-  =  R*  +  <2  ^  2R[ijj  ,i„  ecos*  +  sin  esint  t  cose]  (G.2) 
The  values  of  the  components  of  /  in  units  of  earth  radii  are 


=  0.0577 
=  -0.0324 
=  0.0184 

It  should  bo  notod  that  barred  quanlltios  are  In  units  of  earth  radii  throughout  this  seellon 
6.  2  EQUATORIAL  ORBIT  INJECTION 


mple  eases  are  generated  by  assuming  equatorial  salollllo  orbits  and  eonslderlng  the 
dlsp  need  dipole  to  be  parallel  to  the  spin  asis  and  ,n  the  earth's  equatorial  plane.  The 
problem  then  assumes  the  two  dimensional  geometry  shown  In  Fig.  G-3. 


I  he  satellite  orbit  is  given  by 


R  =  ^  ~ 

1  -  €  cos  ( 4)  -  5  ) 


(6.3) 


where  a  and  e  are  the  elllplle^rblt.semlmajor  axis  and  eeeentrielty 

and  G  ,s  the  major^xls  phase  angle  Introduced  by  the  earth's  rotation' 
becomes,  with  0  =  90° 


respectively, 
Equation  (6.  2) 


2 

r 


R 


+  i 


(6.4) 
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T 


Fig.  6-3  Two-Dimensional  Coordinate  ^stem  for 
Equatorial  Orbit  Calculations 


Differentiating  Eqs.  (6. 3)  and  (6.4),  and  equating  r  and  R  to  first  order,  one  gets 


Ar 


A4> 


+  cos  4> 


+  1. 


c  sin  (4»  -  6) 

1  -  e  cos  (4.  -  6) 


[■ 


'x®'" 


4>  cos 


•]| 


(6.5) 


Several  cases  have  been  considered  for  equatorial  orbit  injection.  These  are  obtained 
from  Eq.  (6.5)  by  applying  special  conditions  which  are  described  below. 


i 

1 


LOCKHEED  AIRCRAFT  CORPORATION 


6-5 


MISSILES  and  SPACE  DIVISION 


Ci!^.  The  0,  orm  »„d  invarlnnl  surface  occur,  at  apogee  when  the  major 

aals  of  the  orb.,  coincides  w.U,  Urn  dlpolo  dlsplacomon.  vector,  which  lies  a.  150. 9-  E 
long.  Then  *=0  5  s,„co  s.n,=,  /2  on  too  avorago. 


A*.  =  2,,  and  o  c  S  iii  a 

rjiirrr'  ™‘““ “““  ■"  "■ 


For  circular  orblls,..  Is  Instructive  to  consider  the  case  of  the  orbit  Inter- 
secllng  the  Invariant  surface  a.  a  point  approx.mate.y  a.  right  angles  to  toe  dipole  d.s- 
P l  e men  vector,  n.s  can  Umn  be  compared  w.to  the  tongen.  case.  Prom  Z  <5. 5, 
e  -  0  .  toe  extreme  values  of  a*  are  obto.ned  by  finding  toe  extremes  of  Its 


coefficient. 


di’  [~^sX  ® 


The  maximum  coefficient  is  at  4,=  60.  O^E,  where 


»?2  =  *  “  60.  9^  .  Then 


A«I-2  =  27/2  and =  0.066y 


(6.7) 


The  tangent  case  Is  again  at  4,  =  l.,o.  9*  E  In  which  case  5  /2  =  0. 0348  „  V2  11.0 

maximum  injection  time  exceeds  the  minimum  by  about  40  percent.  '  ' 

Cas^.  injection  a.  apogee  In  an  elliptic  orbit  Is  sh,r.es.  when  toe  orbit  major  axis  lies 
normal  to  the  dipole  displacement  direction.  For  this  case  5  =  00.9-E  and 
»?3  4  = 'I*  -  60.9*. 
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Then 


2°  '  *  0.056l(y)  +  O.OGGsj 

"  ["^fT^)(T)  *  0-056l(^)  -  O.OGGsj 


(6.8) 


^•*3  =  ’^3  ^  \ 


(6.9) 


Case_4.  Injection  from  the  satellite  at  a  position  other  than  apogee  will,  of  course, 
reduce  the  length  of  the  Injection  time.  This  reduction  Is  obtained  here  for  the  simple 
case  of  an  elliptic  orbit  In  the  equatorial  plane  of  a  geocentric  dipole.  For  R  equal 
to  80  percent  of  apogee,  the  shell  thickness  Is  decreased  by  (0. 8)^  because  of  Increased 
magnetic  field  strength.  The  Injection  position  *  Is  obtained  from  the  orbit  equation 
with  R  »  (0. 8)a  (1  e ) ,  so  that 


0.8 


1  -  € 

1  -  €  COS  $ 


and 


e  -  1  +  €  .  -1 

T  “  1  -  €  (®'®)  sIn$A$^  (6,10) 

6.3  POLAR  ORBIT  INJECTION 

Injection  from  a  polar  orbit  satellite  is  another  simple  case.  Because  of  the  increased 
complexity  introduced  by  the  dipole  field,  the  assumptions  are  made  that  the  dipole  is 
geocentric  and  that  the  rotation  of  the  earth  during  injection  may  be  neglected.  The 
dipole  coordinate  system  is  shown  in  Fig.  6-4,  where  ip  is  the  magnetic  latitude. 
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The  electron  shells  are  thinner  away  from  the  equatorial  plane  because  of  the  Increased 
magnetic  field  strength.  Hence,  the  varlaUon  of  shell  thickness  with  latitude  over  the 
region  of  Intersection  must  be  considered  In  addition  to  the  satellite  position  relative  to 
the  Invariant  surface.  Furthermore,  the  shell  thickness  and  satellite  displacement  are 
now  measured  at  an  angle  »/2  -  y  to  the  radius  vector  r  . 


The  allowable  Increase  In  shell  thickness  normal  to  the  shell  is  given  by 

_±_  =  (1  3sln^;M^^^  ^ 

slny  4  cos^  4 

■  ""orbit  -  ""b 

This  expression  can  be  simplified  by  evaluating  Afp  In  terms  of  which  was  valid 
at  the  geomagnetic  equator.  Bp  Is  a  constant  for  electrons  with  a  given  energy  and 
pitch-angle,  where  p  Is  the  cyclotron  radius.  Then 


Ml  €  U  8ln  i  ^  2  ^ 

(l-ecos^|k)  ®  I 


rr  A 


r 


e 


1 _ 

(1  +  3  sln^i/i)^/^ 


cos  ij) 

n  j.  ,  2  ,.1/2 

(1  +  sin  Ip) 


(6.12) 


Further  simplification  can  be  made  using  the  fact  that  the  magnetic  field  line  and  the 
orbit  Intersect  at  r  so  that 


-  _  a(l  -  _  2 

^  ,  „'i  =  r  cos  Ip 

1  -  c  cos  Ip  e  ^ 


(6.13) 
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These  expressions  allow  Eq. 


(6. 11)  to  be  rewritten  as 


Ar 

4 


3 

n(l  -  c^) 


(6. 14) 


or  rr  ^  ~  ~  -  ^ 

for  mun,  a  ^Tj .T.  “  1“*“  “ 

orbit.  po  ar  or  from  an  equatorial  satellite 


Equation  (6. 10)  may  be  rewritten  as 


sin  « 


(6. 15) 


where  the  ratio  of  injection  radius  to  apogee  radius  is  , 

arWtrarj-  value  of  0. 8  uaed  previously  ta  the  ca  ®"  “  '  "“■or  than  the 

Equation  (8.  M,  ean  be  sImoLled  If  th  .  oonsIderaUon  f  a  i. 

are  roughly  equal,  and  if  r  Is  I  k  '  “‘f  ®q“"torial  orbit  eccentricities 

orbit  apogee.  The  result  Is  ^  'ro®  than,  the  polar 


Ar  .. 

T 


sin  Ip  A^ 

4“ 

COB  4 


Combination  of  the  last  two  equations  yields 


Ml. 


4 

cos  Ip 


(6.16) 


(6.17) 
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6.4  SUMMARY  OF  THE  VALUES  OF  A  4* 


2-44-Gl-l 


Table  6-1  gives  some  numerical  values  of  a*  and  i  for  the  different  eases  dis¬ 
cussed  above.  Satellite  perigee  was  arbitrarily  taken  to  be  at  R  =  1. 05.  and  the 

equatorial  piteh  angle  was  chosen  to  30* .  Several  simple  observations  can  be 
made  about  the  results. 

(1)  For  R  >  5  ,  A*  for  circular  orbits  Is  at  Icastonc  oixler  of  m.'ignitudc 
larger  than  Is  obtainable  at  apogee  from  elliptic  orbits. 

(2)  For  R<  5  ,  A*  for  elreular  orbits  varies  from  one  order  of  magnitude 
larger  to  one  order  of  magnitude  smaller  than  Is  obtainable  at  apogeo 
from  elliptic  orbits. 

(3)  For  elliptic  orbits.  A4*  Is  approximately  two  orders  of  magnltudo  larger 
at  apogee  than  at  0.8  of  apogee.  That  Is.  the  largest  A#  occurs  when 

sin  4>  ^  A4> 

(4)  The  comparison  of  polar  and  equatorial  orbits  Indicates  that  typical 
A4'  away  from  apogee  are  comparable  for  the  two  situations.  It 
should  be  pointed  out  that  one  Important  assumption  Is  made  In  the 
polar  versus  equatorial  orbit  comparison,  namely,  the  earth's  magnetic 
field  Is  not  being  rotated  out  from  under  the  satellite  during  Injection. 

Table  6-1 


SATELUTE  INJECTION  ANGULAR  INTERVALS  AND  ANGULAR  VELOCITIES^^^ 
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6.5  INJECTION  TIME 

To  determine  if  it  is  possible  to  actually  fill  a  given  shell,  the  above  informaUon 
about  A  4.must  bo  combined  with  the  angular  velocity  of  the  satellite  with  respect  to 
the  local  magneUc  field  and  with  the  electron  drift  time.  Table  6-2  lists  some  of  the 
parameters  involved  as  well  as  the  ratio of  injecUon  to  drift  Ume.  11,0  table  is 
based  on  the  orbits  previousl^labelled  as  Case  3.  Figure  6-5  presents  in  graphical 
form  the  dependence  of  ^  on  R. 

Provided  mat  InjecUon  occurs  near  apogee  and  the  angular  velecity  veolors  of  the 
satellite  and  me  earth  are  parallel,  me  behavior  of  may  be  summarized  as  follows 

(1)  For  R  ~  5,  cancellation  of  the  satellite's  rotation  by  that 
of  the  earth  is  the  most  important  factor  determining  the 
size  of/?.  Indeed,  injection  at  apogee  from  an  orbit  of  any 
allowed  eccentricity  will  always  fill  a  shell. 


Table  6-2 

RATIO  OF  THE  INJECTION  TIME  TO  THE  DRIFT  TIME 


Injection  at  Apogee  From  An  Equatorial 
Orbit  Satellite  in  an  Eccentric  Dipole  Field  (Case  3) 


R 

£ 

4>  - 

3  earth 

(milliradians/sec) 

Td 

(sec) 

At.  , 
inj 

(sec) 

^D 

1.5 

0. 176 

0.54 

3,010 

648 

0.22 

2.5 

0.408 

0.17 

1,800 

464 

0.26 

3,75 

0.563 

0.040  ■ 

1,390 

805 

0.58 

5.0 

0.653 

0.0059 

900 

5,050 

5.62 

6.25 

.  0.712 

0.030 

722 

960 

1.33 

7.5 

0.754 

0.043 

600 

737 

1.23 

10.0 

0.810 

0.056 

450 

625 

1.39 
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R  (EARTH  RADII) 
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(2)  Shells  having  R  <  3  cannot  be  filled  fmm  the  relaUvely  high 

oocentrJcity  orbits  ohosen,or  from  circular  orbits.  On  the  basis 

of  the  size  of  ,  elllpUc  orbits  are  about  as  good  as  circular  orbits 
In  the  neighborhood  of  R  a»3. 

■no  ab.v.  co„olu.,o„.  protabl,  ,„r  orbl«  w.U,  repoc.  ,o 

oirlh  8  equatorial  piano  by  anglca  <  0. 5  radinna,  provided  tho  laJcoUon  and 
niaiinoUc  lattitudos  aro  of  tho  aaino  alK„. 

in»J«.lo„  ooours  a,  aomo  toodorato  dlalaaoo  f™„.  apogoo.  tho  valuoa  of/J  will  bo 
much  a,  two  order,  of  magnitude  amallor  than  thoao  appearing  i„  Table  6-2. 

6.0  INFLUENCE  OF  ACCELERATOR  BEAM  FLUCTUATIONS 

In  addition  to  the  radlua  of  gyration  a^  the  thiekneas  of  a  ahell  la  mnuenced  by  the 

apremi  In  aeeelerator  beam  energy.  AE,  and  the  divergence  of  the  beam  Ao  along 

c  I  m,  of  Inleetlon.  ConaequenUy.  the  Initial  thleknea.  Ar  of  a  ahell  at  the  geometrle 
equator  may  be  expressed  as: 


Ar(a  ,  E,  O')  =  2a  + 
o  e 


=  2a  + 
e 


((yAp]  (neCoto’Ao)  j 


=  2a 


((  f 

W  -  mV/ 

+  (cot  Q-Att)  J 

(6.18) 


at  fv  -  ono  .  ''  tnicKness  lor  1-Mev  electrons  injected 

a  erelat  ‘ndlcate  that  ^  and  ^ a  probably 

are  relatively  unimportant  factors  in  determining  Ar.  ^ 
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Table  C-3 


VALUES  OF  -  1  j 


AE 

A  nr 

(Mev) 

(rad 

ians) 

0.05 

0.15 

0.01 

0. 0% 

0.2G 

0.10 

0.11 

0.26 

6.7  MAGNITUDE  OF  INJECTED  FLUX 


An  estimate  of  the  magnitude  of  the  injected  eleetron  flux  F  is  most  easily  obtained  by 
considering  the  number  of  electrons  crossing  tlic  equatorial  plane  per  em  per  see.  The 
radial  thickness  of  the  shell  at  the  equator  is  assumed  to  be  2.5a^;  the  aeeelcralov  injects 
a  current  I  for  a  time  Atj^j  at  a  geomagnetic  radius  Electrons  arc  injeeted  at  the 
geomagnetic  equator  with  a  pitch  angle  a  with  respect  to  the  local  field,  and  thereafter 
cross  the  equatorial  plane  with  the  same  pitch  angle. 


/  total  number  of  \  /number  of  times  per  see  each  eleetron' 

p  _  \eleetrons  injeeted/  \ _ erosses  Uie  equator _ , 

equatorial  /eross  seetional  area  in  equatorial  plane  \ 

plane  \  as  seen  by  the  eleetrons  / 


1  Q 

6.3  ^  10  At,  , 
in) 

27rr^(2. 5)a^  eos  a 


_l _ 

bounee 


1  8 

6.3  ^  10  I 

27rr  (2.5)a  eos  at 
e  e 


-(J 


A'l’ 


satellite 


4> 


earth 


) 


ounee 


(6.  19) 


This  equation  assumes  that  no  injeeted  eleetrons  are  lost,  and  the  eleetrons  are  spread 
uniformly  around  the  shell. 
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Value,  of  F  and  an  eaUmale  of  the  Van  Allen  elootron  flux  above  0. 9  Mev  are  given  in 
Table  G-4  for  tbo  inJecUon  parameter,  of  Caae  3  found  from  Table  6-2. 

Table  6-4 

ALTITUDE  DEPENDENCE  OF  INJECTED  FLUX  AND  ESTIMATFn 
VALUES  OF  VAN  ALLEN  ELECTRONS  BETvJ^EEN  0  S  S^l  MEV 

(Flux  units  are  electrons  |)cr  em^  per  see) 


Van  Allen 


2  X  10 
4  X  10 


The  eorreeted  eleetron  speetrum  of  Walt,  et  al.  (Ref.  17).  was  used  to  estimate  the  Van 
Allen  flux.  Extrapolation  of  their  spectrum  indicates  a  flux  of  eleetrons  between  0. 9  and 
1. 1  Mev  of  the  order  of  10  per  em^  per  see.  Values  of  Urn  eleetron  flux  at  different 
geomagnetie  distanees  ean  be  arrived  at  by  sealing  the  Walt  flux  of  5  x  lo®  eleetrons 
greater  than  20  kev  per  em  per  see  by  the  variation  of  quiet-day  flux  with  distanee  as 

given  by  Dessler  (Ref.  3).  The  Van  Allen  eleetron  flux  so  obtained  is  uneertaln  by 
one  or  two  orders  of  magnitude. 

It  should  be  noted  that  the  IntUal  angular  diatributlon  of  Injeeted  eleetron,  will  be 

different  from  the  angular  distribution  of  Van  Allen  eleetrons.  This  dlreeUonal  feature 

of  the  Injeeted  flux  may  be  useful  In  diserimlnating  against  the  flux  of  naturally  oceurrlng 
eleetrons. 
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Section  7 

BOUNCE-TIME  MEASUREMENT 


7.1  INTRODUCTION 

TIic  major  obstacle  to  the  bounce-time  munsurement  is  tlic  detection  of  tiie  injected 
electrons  after  their  first  bounce.  We  have  proposed  Uiat  tlicse  electrons  be  detected 
at  the  satellite  by  the  radio  signal  from  Uie  electron's  synchrotron  radinlion  near  the 
mirror  point.  To  demonstrate  the  feasibility  of  this  metliod,  \vc  will  compute,  in  a 
typical  case,  the  signal  strengUi  from  tlic  electrons  and  tlic  sources  of  background 
which  would  interfere  with  the  reception  of  the  signal. 

7.2  SATELLITE  AND  ELECTRON  MOTION 

We  will  consider  that  the  satellite  is  at  a  low  altitude  (a  few  hundred  kmsjand  is  in  a 
circular,  polar  orbit.  The  pulse  of  one-Mev  electrons  (100  amps  for  <1  x  10”^sec)  is 
injected  nearly  perpendicular  to  the  direction  of  the  local  magnetic  field.  Aftei  the 
time  for  one  bounce  has  elapsed,  the  electrons  return  and  mirror  at  a  point  near  the 
injection  point,  but  at  a  lower  altitude  and  east  in  longitude.  During  the  bounce  time, 
the  satellite  has  also  moved  in  latitude. 

We  will  restrict  ourselves  to  one  special  case  throughout  the  discussion  of  the  bounce¬ 
time  measurement;  namely,  injection  on  the  field  line  which  intersects  tlie  earth  at  a 
magnetic  latitude  of  45”  and  crosses  the  magnetic  equator  at  2  earth  radii  from  the 
earth's  dipole  center  (r^  =  2  R^.)  .  At  this  position, 

B  =4.6  X  10  ^webers/m^ 

p 

oj  =  2.9  X  10  rad/sec 

r  =  440  Kc/sec 
o 
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=  100  m 
Tjj  =  0. 2  see 

Satellite  veloeily  =  8  km/see 
Longitudinal  drift  veloeify  of 
the  mirror  point  =  13.5  km/see 


After  the  first  bounce, 
shown  in  Fig.  7-1. 


the  relative  positions  of  the  satellite  and  the  electrons  are 


."^\T.TLUV!: 

IG  ICr.i 

iriv'OTiOtf  - -  V 

0  electro:  13 

Fie.  7-1  RelaUve  Position  of  the  Satellite  and  The 

Electrons  After  One  Bounce 


^e  eleetrens  aetually  ndrror  at  a  lower  altitude  U.„n  U,e  Injeetlon  point  for  any 
jeetion  plteh-angle  o,  other  than  90-.  Of  course.  If  o  were  exactly  90- ,  the 
inteeted  electrons  would  strike  tl.e  satellite  after  one  eyelot,-on  revolution  at  the  time 
oflnjeefon.  However,  If  89- this  will  „o,  oeeur,  and  we  will  assume  that 

w  s  etween  80  and  89-.  An  expression  can  be  derived  for  the  difference  in  altituc 

the  cenT'’ no  “i  '  f 

center  of  the  earth  to  the  injection  point  be  r  ,  and  to  the  mirror  point  r  Th. 

m  " 


r. 

1 


Rp  +  h. 
E  1 


(7. 
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where  Rj,  is  Ihe  radius  of  ihe  earth,  ami  h.  and  h^^^  are  the  alliliides  of  the  injection 
point  and  mirror  |K)inl,  respectively.  From  Ujc  inatjnelic  momeiU  invariance,  vve  have 


Bi  ==  B,^^  sin-  o-j 


(V.2) 


From  Eqs.  (2.  i)  and  (2..*}a),\vc  fiiul  Unit 


If  wo  take  case  hj  ,  h^  <<  ,  then 


(7.3) 


(7..I) 


where  c  ^  is  defined  by  =  |  -  c.  .  £^<10“.  Tl,c  result  does  not  depend  on  Ute 

value  of  r  to  the  first  order. 

6 

7,3  SYNCHROTRON  RADIATION  FROM  ELECTRONS 


The  instantaneous  power  radiated  by  an  electron  witli  velocity  v  and  total  energy  E  , 
that  is,  kinetic  energy  plus  tlie  rest  mass  energy,  is  (Ref.  21): 


(7.5) 
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In  Ui.  canc  of  an  election  apirallng  alon*  a  field  line  with  pitch  angle  «  , 


2  3 

a  o  p  (a 

P  =  (6  X  10®)  - - £.  sin^ 


•(ij 


(7.6) 


where  p  .  ,/c.  „d  “«  ‘o  "•  •’■n  cyclotron  Migular  frequoncy  and  qrclotron 
r^e,  rea^cuvcly.  The  power  radlaled  can  then  he  expreaeed  In  lerm.  o7th. 
local  magneUo  field  strength  and  the  piU;h  angle. 


P  =  (6  X  10®) 


& 


c®  B®  sin® 


(7.7) 


If  P-  1  ,  then 


P  »  (1.6  X  10"*^)  B®  sin® 


(7-8) 


Onc-Mev  electrons  mirroring  near  the  surface  of  the  earth  at  a 
45* ,  results  in 


magnetic  latitude  of 


P  =  3  X  10  ®®  watts 


In  the  pulse  of  electrons  there  are  2  ‘i  x  alone..  . 

tho  ai.i  j  u  -6  *  electrons,  so  that  the  power  radiated  bv 

the  pulse  Is  about  10  watts  at  the  mirror  point. 

Zr  “‘“•■■‘'’"‘'cn  cd  an  angular  distribution  which  Is  not 

iTcc  l  r‘"T 

n  accelerating  electron,  and  has  applied  the  general  formula  to  the  case  of  a  high 
gy  electron  (fl  -  1,  moving  with  constant  velocity  In  a  circular  path.  Since  the  motion 
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Is  U.0  . . . os  . . . .  (o.vololron)  f,oq.,o„ov 

^  ■  ko/soc.  I'-tnllnstaMncoiis  ra,||,,,,,,l  |„io  IIK!  I,'''  |,,,,.,,,o„,o  |„  ' 

n«U  orlllcol  harmonic  number  In  .lofinoci  I,,  =  3/2  (E/mo^  f  .  oaao .  "  ’ 

=  40  .  If  fho  condition  n  «n^  is  sniisfled.  then 


P„  =  (4  X  10®)  c^/c  o)^ 


If  llio  condition  n  »  1  is  satisfied,  then 


=  (3.72  X  10®)  c® 


•  9  /  c-  \2  r” 

’/c  u  — \  JL  I  K  / 

c  (  2  1  2  J  ,  r)/3 

\mc  )  II 
\  /  c  c 


(•/)  d»/ 


(7.10) 


where  Inlhe  Alryinlegral.  Tmnlx.nll»n  and  llnrlmnn  (Ref.  22)  have  evnlualml 

numerleaU,  .he  Integral  In  Mq.  ,2.  P,,„oe  2-2  aho...,  the  values  ef  P  „„ 
of  frequeney;  .he  flra.  seven  harmonies  are  given  oaplleltly,  bu.  for  higher  n's,  only 
re,.resenu.tlve  values  of  P__  are  given.  The  most  Interesting  feature  of  speetrum  Is 
he  fact  that  due  .0  relativis.le  effeeb,  the  maximum  ,x,„er  Is  radiated  In  al»n.  the  2o"' 

11“.’  ““'W 

Al»ve  2  me/see.  U.e  density  of  the  lines  of  the  dlserete  spectrum  becomes  sufficiently 

ralteraT  ”  ^  •<"»  "-e  power 

..adlated  at  frequeney  a  In  the  Interval  dr.  wc  find  the  eentinuous  speetrum  to  be 


P(t')  =  P 

"c  " 


(7.11) 


This  spectr.il  distribution  is  shown  in  Fig.  7-3. 
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■itc.  1-2 


ting  w!"h  a'Frequwicy  0^40X76"^  “"‘ThjTfi  Eleclron  Circula- 

expncity;  at  Higher  freqaLl^ia 
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FftEOUEJJCY 

Fig.  7-3  The  Continuous  .Appro.ximation  to  the  Spectrum  of  Synchrotron  Radla 
tion  from  an  Electron  Circulating  with  a  Frequency  of  440  Kc/sec 


2-4‘l-61-l 


The  angular  distribution  of  tlie  radiatcci  |)o\ver  has  also  been  derived  by  Schwinger. 

Tliu  expression  is  complicated,  and  for  the  purpose  of  this  study,  it  will  be  sufficient  to 
use  the  rough  characteristics  of  the  distribution.  Namely,  the  power  is  a  maximum  in 
the  orbital  plane  of  the  electrons  and  falls  off  rapidly  with  increasing  values  of  the  angle 
defined  in  Fig.  7-4. 


B 


Fig.  7-4  Angular  Coordinate  for  Synchrotron  Radiation 


Virtually  all  the  radiation  is  emitted  at  angles  less  than  =  (l/jr  .  In 

our  case,  y^^^  ^  20“  .  If  we  assume  that  the  angular  distribution  of  the  radiated  power 
is  constant  for  0  <  y  <  ,  and  zero  for  y  >  y^^^  ,  then  we  can  compute  the  power 

incident  on  an  antenna  at  a  distance  L  from  the  circulating  electrons,  where  L  is 
much  larger  than  the  electron's  cyclotron  radius. 
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2  ii-r.1-1 


8 


N{i)Ar 


NAP(»)t\r 

*» 

•Ib  L~  y 


(7.12) 


S(i')Ai'  Is  the  sluiKil  Hlronf^th  at  Uie  Input  to  a  nullo  receiver  which  Is  lune<l  to  a  rrc<|iioncy 
V  aiKl  has  a  liaiuhvllh  Ai*  i  A  Is  the  effective  area  cif  the  antennae  and  N  ilie  number  of 
clrcuhitluK  clcelrons.  If  we  lake  the  followliiK  parameters  for  the  satellite  ex|M*rlmcnl. 


A  =  1 

e  =  l.T)  Mc/sec 
At'  =  1  Mc/soc 

P(i')  =  13  X  10  wntls/Mc/sec 

N  =  2.. '5  X  lO*'*  electrons 

y  =  0.3  r.adians 
max 

b  =  3  km 

-If* 

then  the  signal  strength  S  Is  10  '  watts. 

7.  I  SYNCHROTRON  RAOIATION.SIGNAL  DURATION 


Before  examining  the  soiirees  of  baekgrouiul  noise  whieh  effect  the  reeeption  of  the 
syiiehrotron  radiation  signal,  an  estimate  will  have  to  Ikj  made  of  the  time  duration  of 
the  signal.  On  one  hand,  we  would  like  the  signal  to  be  as  short  as  {wssible  to  minimizf! 
uneertalnty  In  the  measured  bouiiee  time;  on  the  other  hand,  in  so  doing,  the  masking 
effeet  of  the  background  noise  becomes  prohibitive.  A  calculation  will  show,  however, 
that  the  duration  of  the  signal  is  relatively  insensitive  to  any  practical  control  we  might 
attempt  to  impose.  If  the  satellite  has  a  directional  antenna,  such  as  a  parabola,  the 
signal  duration  is  determined  by  the  length  of  time  required  for  the  eleetrons  to  pass 
through  the  acecptancc  angle  of  the  antenna.  The  relative  position  of  the  satellite  and 
the  eleetrons  ai’e  shown  schematically  in  Fig.  7-5. 
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a 


Fig.  7-5  RelaUve  PosiUon  of  the  Satellite  and  Electrons 


Since  the  electrons  are  radiating  direcUonally.  there  is  the  quesUon  of  whether 
e  8  g8<d  duration  would  be  dotermlaed  bjr  the  bearawidth  ol  the  antenna  or  by  the 

beamwldth  bait-angle  to  be  equal  to  ;  the  signal  duration  will  Uien  be  given  by 


(Rj.ef/3 


Ly, 


tl/2 


max 


)  -(Hr 


cf/S  -  Ly 
^  max/ 


(7.13 


ir/Z  -  €j  .  Equation  (7.13)  assumes  that  the  electrons 

r  7.3  r  “  "’•»  "»•  then 

r<q.  (7. 13)  becomes 

vl/2 


,ef/3 


+  L  y 


max 


jl/2 


(7.14) 
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ll-Oi-l 


'l  lii*  l•\|Ir(^s.s^(ms  for  lake  Into  iK'coiiiit  iho  fact  llial  ••Icclrons  {vis;;  tho  :iiitoiiii:i 
l\v\vv  as  tlioy  mlrmr  and  Iwimco.  Tin-  sljtiial  thirnllon  is  nlvcn  In  llio  l(*ll..vviii{r  lahlo 
for  a  few  liijoolinn  iiUrIimIukIos. 


fi  X  10  *  sec 
1.5  X  10  sec 
8  X  lo"'*  sec 

It  Is  not  pos.sihic  to  Inject  with  €  j’s  much  Kfcaler  than  lO”  without  lowering  the  mirror 
IkjIiU  to  such  a  degret;  that  serious  scattering  losses  occur  at  the  first  bounce  In  the 
other  hemisphere  licfore  the  electnais  return  to  he  monilcn-ed.  ■Iliere  will  he  an  Intrinsii: 
beam  spread  determined  by  the  electron  accelerator  Itself,  so  that  It  will  not  be  po.ssibic 
to  define  tluj  injection  pltcli-anglc  to  licltcr  than  a  few  degrees,  ^•ol•  the.se  reasons,  the 
duration  of  the  .synchroton  r.adlalion  signal  will  be  about  10"‘'  see  In  most  e.ascs. 

The  lx)uncc  time  in  our  special  case  is  0.2  see.  The  experimenl.’il  uncertainty,  there¬ 
fore,  associated  with  the  lx)uncc-timc  measuremcait  would  lx*  alwiit  0. 1  percent. 

7.5  BACKGllOUND  NOISK  SOURCES 

There  are  a  luimhcr  of  sources  of  radio  noise  which  radiate  in  the  10  Me/scc  region, 
and  will  ix!  picked  up  by  the  receiver  on  the  satellite.  ‘J’lie  ai)ility  to  detect  the  injected 
olectroii.s  will  be  limited  only  by  the  noise  background;  and  the  signal-tn-noise  ratio  must 
be  estimated.  We  will  consider  the  following  sources  of  noise: 

o  Thermal  emission  from  the  galaxy 
o  Thermal  emission  from  the  sun 
o  Radiation  from  the  V'an  Allen  electrons 
o  Intrinsic  receiver  noise 


0 

5“ 

10* 
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.oU,or  possible  source  of  baekKround  noise  Is  radio  U-ai.sralsslo„  from  Uio  carUi.  This 
«11l  presumably  bo  negligible,  however,  since  the  Ionosphere  Is  o,.a,p,e  to  radio  fre- 

.luoneles  less  U.an  about  30  Me/see.  Most  of  .he  following  h, formation  about  radio  no.se 
lias  been  taken  from  Pawsey  and  Brneewoll  (Ref.  23). 

7.5.1  Gnlaetie  Noise 


Ihc  radio  omission  from  the  gniaxy  can  bo  representetl  ns  coming  from  a  blaekbody 
with  temperature  T|^ .  Radio  astronor-ers  scan  the  galaxy  and  present  tlieir  results 
ns  a  eolesllal  mapping  of  T,, .  From  the  blnekbotly  radiation  law.  the  ,x,wer  from  a 
galactic  distributed  source  striking  an  antenna  which  has  a  bcamwidlh  talf-anglc  of 


^inax  •’cpresenled  by 


Pbg(«')Ai;  =  (9.67  X  t  Ayl_  ai» 


nin.\ 


(7.15) 


whore  A  Is  Uic  area  of  the  antenna,  and  Uie  satellite  radio  receiver  Is  tuned  to  a  fro- 
queney  v  ai\<i  has  a  bandwidth  Ai'  (boUi  in  Me/see). 

Wlien 


Then 


V 

Ar 


A 

y 

max 


10  Me/sec 
1  Me/see 

1 

20“ 


1. 1  X  10~^®  T. 


(7.16) 


Pawsey  and  Braeewell  give  the  values  of  Tj^  as  a  function  of 
greater  than  20  Me/sec.  Extrapolation  of  their  numbers  lead  to 


V  for  frequencies 
galactic  temperatures 


7-12 


LOCKHEED  AIRCRAFT  CORPORATION 


MISSILES  and  SPACE  DIVISION 


al  10  Mc/soc  of  about  10^*  K  In  the  bright  center  of  the  galaxy,  and  lO**"  K  In  the  cold 
parts  of  the  galaxy.  The  inaxtinum  galactic  noise  Is,  therefore, 

P|^  =  10  watts 


7.r>.2  .Solar  Noise 


If  Eq.  (7. 15)  is  multiplied  by  the  ratio  of  the  solid  angle  snbtendcti  by  the  sun  at  the 
|)olnl  of  observation,  to  the  solid  angle  defined  by  the  antenna  bcnnuvidth,  and  if  Tj^ 
is  the  effeetive  solar  radio  tem|>erature.  then  the  solar  noise  |)ower  is  represented  by 


(1.1  X  10 


(7.17) 


n  -  7  X  10  steradians  and  T,  =  10***  K 
s  b 


Therefore 


P,  =2x10^®  watts 
bs 

7.5.3  Noise  From  Van  Allen  Eleetrons 

The  trapped  Van  Allen  electrons  will  emit  synchrotron  radiation,  and  the  spectrum 
of  this  radiation  depends  sensitively  on  the  Van  Ailen  ciectron  spectrum,  particularly 
in  the  energy  interval  around  one  Mcv.  Since  the  ciectron  spectrum  is  not  known  wcil, 
a  calculation  of  the  synchrotron  noise  would  not  be  possible.  However,  Kraus  and 
Ko  (Hcf.  24)  have  scanned  the  region  of  the  Van  Allen  belt  with  a  ground  station  radio 
receiver.  They  report  that  the  effective  blacklx>dy  temperature  in  this  region  is  about 

5 

10  ®K  (extrapolated  to  a  frequency  of  10  Mc/sec).  Applying  Eq.  (7. 15)  one  finds  a 
value  of  Van  Allen  electron  synchrotron  radiation  noise  to  be  10“^^  watts. 
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7.5.4  Receiver  Noise 


“s'  ““T -»■«<- 

nrs.  su.«e  or  the  .1,1“  »«> 

..olsc-frco  soever  »„o,o  «.o  si^!!!  .Trolls"'’"""'^  “  "  ""  " 

P,,  =  0.37x.„-n,^^__^_  ,r.,„ 

of  N  =  1  for  Uie  best  nmaon#  su  .  ^  ”  ''•'‘lue 

uic  oesi  present-dsiy  receivers.  For  t  =  30o»  k 

^br  ~  ^  ^  watts 
7.6  SIGNAL- TO- NOISE  RATIO 

The  dosorlpllon  of  Ihc  background  . .  a  , 

galactic  noise  which  frivno  n  «  «  *  *  ^  dominant  source  is 

noise  Which  gives  a  noise  signal  strength  at  the  saiPllifo  ,af  ,n“l4 
Is  ten  times  the  strenirfh  nf  »h  i  ^  satellite  of  lo  watts,  which 

es  me  strength  of  the  signal  we  are  attempting  to  cleteet  Thi=  i 

howovor,  that  ,1,0  slgnal-to-nolso  ratio  Is  1/10  *  k.  !  , 

value  without  affootlng  the  rooootlon  of  111  d  ,'  t  """ 

constant  and  can  bo  biased  out  Onlv  fl  “  ™  '“"6  “"C  taokground  Is 

--  with  the  signal  dm  .rirrirn  “ 

random,  then  over  n  limp  t  n  •  -we  assume  the  noise  to  be 

average  value  of  deviation  from  its 


(rAr) 


1/2 


U.iy; 


:iz‘.  rr;r.:.r:“':  “-.rtr 

noise  ratio  of  one  which  hop  *  k,  •  ,  ’  ™  a  signai-to- 

>  IS  acceptable  in  determining  thn  Qimr.if> 
the  electron's  presence.  ^  simple  yes-no  situation  of 
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7.7  BOUNCE- TIME  EXPERIMENTS 

In  most  cases,  the  results  of  a  bounce-tIme  measurement  cannot  bo  estimated  before¬ 
hand  because  of  Uio  lack  of  knowledge  about  tlio  geomagnetic  field  o.'ccopt  iu  very  idealized 
models.  One  ealculation  has,  imwever.  been  made  of  tl>e  bounce  time  predicted  by  two 
different  g<K):nagnelie  models.  The  purpose  of  the  eomparlson  is  to  determine  whether 
the  experimental  aecuraey,  which  will  about  0.1  percent, is  sufficient  to  easily  dif¬ 
ferentiate  between  the  models.  The  two  models  we  have  used  are  the  dipole  model  and 
Welch  s  512-torm  harmonic  expansion  (Ref.  2G)*  The  results  of  the  comparison  are 
shown  in  Table  7-1  for  four  different  Injeetion  points.  One  point  is  near  tl>e  Capo  Town 
anomaly  where  tli(  .:iffcrcnoo  in  pre<liotcd  bounoo  times  should  be  nearly  a  maximum. 

The  difforcnoc.s  ai'c  about  10  percent,  and  so  the  o.xpcrimontal  aoouraoy  is  quite 
adequate. 
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Section  8 

CONJUGATE-POINT  MEASUREMENT 


(U-l 


8.1  INTRODUCTION 

The  critical  question  in  tiic  conjugate-point  incasurenicnt  is  the  fonKlbility  of  detecliiig 
the  electrons'  presence  in  the  upper  ntniosphcrc  by  rndnr  reflection  from  tiie  trail  of 
ionization  produced  by  Uie  eleetroiis  as  Umy  splrtU  down  a  field  line.  In  a  model  which 
will  prove  to  be  too  simple  to  describe  Uie  real  sitiKition,  Uie  electrons  are  confined  to 
a  narrow  pencil  with  a  cross  section  of  a  few  squai'e  centimeters  defined  by  tlic  accel¬ 
erator  at  the  time  of  injection.  This  narrow  beam  will  enter  tlie  upper  atmosphere 
normally,  will  ionize  but  not  scatter,  and  will  stop  at  Uie  end  of  its  range  at  an  altitude 
of  about  GO  km.  The  atmospheric  (iensity  increases  exponentially  as  Uie  electrons 
penetrate  to  lower  altitudes;  also,  the  energy  loss  per  collision  increases  as  the  elee- 
trons  slow  down.  Tlierefore,  virtually  all  of  Uie  atmospheric  ionization  occurs  near 
the  end  of  the  electrons'  range,  and  wo  will  assume  Uiat  the  ionized  trail  is  confined  to 
an  altitude  interval  between  GO  km  and  70  km.  Since  the  1-Mov  electron  accelerator 
injects  2.. 5  x  10  electrons  in  one  pulse,  Uie  energy  contained  in  Uic  beam  is 
2.Ti  X  10  ev.  It  requires  about  50  cv  to  pi'oduce  an  ion  pair  in  a  collision,  which 
means  that  over  the  lO-km  ionized  trail  Uicre  arc  5  10^^  ion  pairs,  or  an  ion  p,air 

linear  density  of  5  x  10  ion  pairs/ni, 

Pawsey  and  Brascwell  (Ref,  23)  discuss  Uie  problem  of  radar  reflection  from  Uic  ionized 
trail  of  meteors  traversing  Uic  upper  atmosphere.  The  eharacteristics  of  Uie  meteor 
trail  are  very  similar  to  those  of  the  electron  trail,  namely,  the  ionization  is  eonfined  to 
a  peneil  of  a  few  square  eentimeters  in  eross  section  and  occurs  mainly  at  altitudes 
around  80  km.  A  bright  visual  meteor  will  produce  about  10^^  ion  pairs  per  meter  of 
path,  and  is  easily  detected  by  radar  reflection.  Tlie  electron  beam  produces  a  trail 
with  an  ion  density  an  order  of  magnitude  larger,  and  sliould  therefore  give  a  substantial 
radar  reflection.  Because  of  the  effects  which  arise  in  ionized  columns  eontainiiig  very 
great  ion  pair  densities,  this  conclusion  is  not  necessarily  valid.  The  reader  is  referred  to 
the  chapter  on  meteor  detection  in  Pawsey  and  Braeewell. 
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In  11., s  Sh,<iy  wo  i,avo  rosO  lolci  oursolvos  ,o  ,l,o  dolooUon  of  U,o  electron  trail  by  radar 
r.-neeliea  melboda.  However,  addltlcntl  enleulalioas  ...ay  show  that  U,o  trail  is  opUeally 
visible  as  well.  It  is  not  possible  to  eoiKludc  immedlaloly  that  sineo  a  meteor  trail  Is 
visible,  Ihe  eleetron  fall  wlU,  a  comparable  Ion  pair  density  will  be  visible.  In  the  ease 

or  a  moleor,  the  amount  of  lls  cnerRy  which  rocs  Into  light  Is  thought  lo  be  determined  by 
«  different  mechnnisni  Ilian  in  the  electron  ease. 

The  above  model  of  Un-  electroa  Ionised  tn.il  Is  Inadequate  because  U,e  oleetrons  are 
not  eonflned  to  a  narrow  pen.  il,  and  llm  effect  of  sealiering  as  U.e  electrons  slow  down 
In  Uie  atmosphere  Is  npprcclnhle.  inUie  follovvIiiK  sections  we  will: 

(1)  ConipiHo  the  shape  of  the  electron  Iienm  as  It  strikes  the  top  of  the  atmosphere 

(2)  Determine  the  effect  «»f  senttcrlii}; 

(2)  Compute  the  radar  reflection  from  an  c.xtended  rcRlon  of  Ionization 


BEAM  AT  THE  CONJUGATE  TO  THE 


INJECTION 


The  shape  of  the  eleetron  beam  after  one  tr.aversal  from  the  Injection  point  lo  the  conjugate 
point  will  be  determined  by  several  effects. 

8. 2. 1  Magnetic  Moment  Nonconservation 

Elect, -oas  which  .are  Injected  .at  ve.-y  small  pitch  angles  In  a  curved  nnagaetie  field  do 
not  satisfy  U,e  ooaditlons  for  .adhabatle  Invarlanee  of  the  magnetic  moment  so  that  their 
pilch  angle  .at  Ub  point  conjugate  to  the  Injection  point  Is  not  equal  to  the  Injection 
pilch  angle.  This  is  because  at  small  pitch  .angles  appreciable  eh.angc  la  the  field  ean 
occur  within  the  distance  covered  by  the  guiding  center  during  one  eyelolron  revolu¬ 
tion.  The  effect  of  the  field  curvaltiro  Is  such  that  electrons  which  arc  Injected 
p.arallcl  to  a  field  line  will  arrive  at  the  conjugate  point  wiU.  a  finite  pilel,  angle  The 
following  is  a  correct  descrtpllon  of  Uic  motion.  An  electron  would  move  parallel 
to  the  field  line  until  the  field's  curvature  causes  7  and  B  to  diverge,  which  creates 
a  non-zero  plteh  angle.  Tl,e  eleetron  would  U,en  have  a  component  of  velocity 
perpendicular  to  U,c  field  and  would  start  to  spiral,  keeping  the  center  of  gyration 
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cn  the  orlplinl  Injeeliou  liel.l  ||m..  n,..  n,.,.!..  .  ..Mih„.e  |„  inon  .iM..  .n,ii|  (he 

|iolnl  is  rebelled  where  (he  iiiiiiiteiil  liiv.'ii  l;iiiei>  is  ulieyed,  :incl.  (Iiei'e.irter,  the 

inetlon  wlJI  he  properly  .leseril.ed  hv  lh:.l  low.  AUviUt  (Itef.  {,  )  ,||scussos  this  sllii.itlon, 
nnd  sets  (he  followluj;  erllerjoH  f«r  (he  iioiKncufe  inoiiieiii  e<iiist!rv:illon.  The  perpeiidlenhir 

eoi»|«,no»(  of  (lie  eh?elroii'K  veloeMy  l,e  nnieh  m  .  i,!,.,-  (h.in  (he  olecti  fni's  loiiKl(ii(lln;il 
drift  velocity. 


V  ::iii  O'  >>  v,. 

•>  (8,1) 

This  condlllon  on  nr  liniKises  the  KicrKest  re.striclloii  :i(  (he  e<pi;i(or.  If  we  orliMriirlly  set 
Vjj  =  (0.01)  V  sin  ,  ;ind  If  «  (.  (hen 


•Ip  .  2 

^  ell^- 
c?  e 


(8.2) 


(•or  (he  field  line  which  erossos  the  m;i|;ne(ie  ef|ii.’ili»r  .'It  2  e.'irlh  rodil,  o  =  l..'»".  An 

electron  will.  i.  pitch  n.mie  nl  (he  cpinlor  will  . . .  n  field  line  to  the  top  of  the  ntmo.s- 


phcrc*  where  its  pilch  nni^le  will  he  ir  . 

o 


•»  '» 

K  fr 

M  li  (»••'•) 

o  V 


For  the  field  line  r  =  2  earth  radii  and  o  -  1.5",  ihuii 
e  e 


(V  ~  5 
o 


The  eyeh.lron  radius  at  the  top  el  (he  almesphere  will,  therefore,  be  about  lOni. 


i) 
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8.2.2  Aeeolcrator  Beam  Spread 


Allhougl,  Ihc  I, cm  Is  well  colllmalcd  by  U>e  elcclron  eccoleialor,  U>oro  will  be  some 
rinlle  beam  .llveigenee,  Presaraably,  u>ls  beam  spread  will  be  less  Umn  5-  The 
eoneluslon  whieh  ean  be  drawn  from  U.e  arpimenls  In  the  abovo  paragrapl.  Is  11, at 
an  olcclron  injeoled  with  any  pilch  nnRlc  less  Uian  5*  will  arrive  at  tlic  conjugate 
|x.int  with  a  pitch  angle  of  about  6\  Therefore,  U.e  effect  of  U,e  beam  spread  at 
Injcetion  is  completely  masked  by  Uie  magnetic  moment  non-conservation. 

8. 2.  .1  Energy  Dependence  of  the  Longitudinal  Drift  Velocity, 

The  injected  electrons  will  have  a  finite  energy  spread  determined  by  the  accelerator 
(3  imreent  to  ■!  percent  for  U,e  Varlan  Ihtlse  Transformer).  The  longitudinal  drift 
voloelty  of  U.e  electrons  is  energy  dependent  ConsequenUy,  U,c  electrons  will  arrive 
at  the  conjugate  point  wiUi  a  longitudinal  spread.  EquaUons(3. 13a)  .and  (3.14)  give 
the  rate  of  change  of  the  electron's  longitude  averaged  over  one  bounce. 


d</>  3  p  V 

dt  “  2  2  „  (8.4) 

e  r  B 
e  e 


An  electron  with  momentum  p  and  velocity  v  shifts  in  longitude  an  amount  A  4.over 
U.e  tln.e  ,-cquired  to  move  from  tl.e  injection  point  to  the  conjugate  point. 


A  0 


6 


_ B_ 

er  B 
e  e 


(8.5) 
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At  the  conjugate  point,  this  shift  in  longitude  would  correspond  to  a  linear  displnccinonl 
given  by 


f  • 


0 


o  B 

e 


(8.r>) 


where  r^  is  in  units  of  earth  radii.  We  want  to  compute  the  change  in  resulting 
from  a  change  in  the  electron's  kinetic  energy  T. 


L  _  6  /_  \-3A  1 

~  rr  (•'e)  2  -  (8.7) 

For  a  1-Mev  electron  moving  along  the  field  line,  r  =  2  earth  radii,  and  for 
dT  =  (0.03)  T  =  3  X  lo'^ev,  ** 


dDj^  =  50m. 

8.2.4.  Beam  Blow-up  Due  to  Coulomb  Repulsion 

Consider  a  beam  of  electrons  of  length  L  with  cross-sectional  radius  b  moving  with 
velocity  V  and  containing  N  electrons  distributed  uniformly  through  the  beam. 

If  b  «  L,  then  at  the  surface  of  the  beam,  in  a  region  not  too  close  to  the  ends,  there  is 
an  outward  directed  electric  field  given  by 


r 

! 

1 


e 

bL 


(8.8) 


I 

1 


LOCKHEED  AIRCRAFT  CORPORATION 


8-5 


MISSILES  and  SPACE  DIVISION 


2-44 -G 1-1 


Duo  .0  tho  magnetic  field  aaeoclated  «1U.  Ihe  moUon  of  the  eleotrone, 
inward  directed  force  at  Uio  surface  given  by 


Uiere  is  aJso  an 


vVn 

bL 


The  ratio  of  these  opposing  forces  is 


(8.9) 


or 


«« V 

o 


(8. 10) 


(8.11) 


Since  v  can  never  be  identically  equal  to  c  tho  r  , 
otolme  the  beam  would  blow  up  m  our  ca.I  ,  ! 

be  th<^ht  of  as  a  highly  polarlsable  medium.  The  beam  would  blow  up  until  a  point  Is 

m  il  ?  ‘nX'vidual  eTectrl 

distan  hi  I  ““  ^'P“'olon.  An  argument  can  be  made  as  to  the 

stance  which  must  exist  between  electrons  before  they  essentially  lose  communication. 

following  sketch  (Pig.  8-1)  represents  the  situation  where  two  of  the  beam  electrons 

are  separated  and  neutralized  by  the  polarized,  fully-ionized  plasma  between  them 
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PLASMA  CHARGE 
DENSITY 


Fig.  8-1  Separation  and  Neutralization  of  Two  Beam  Electrons 
by  Polarized,  Fully  Ionized  Plasma 


In  the  interval  h  there  is  a  net 

be  assumed  that  n  is  uniform. 
6 


negative  charge  density  n^,  and  for  simplicity, it  will 
In  this  region.Poisson's  equation  for  the  potential  is 


(8.12) 


Integrating  twice,  we  have 


(8.13) 
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In  onlor  for  the  pinsmn  to  hove  liolorteed  to  tho  oxicnt  shown  in  the  sketch,  the  electric 

force,  .mst  have  everce.no  the  thermal  moUen  of  the  plasma.  If  thi.  I.  true,  then  the 

^  n  a  in  Eq.  (8. 13)  muat  correspond  to  a  change  In  energy  of  at  least  1/2  kT,  where 

is  the  temperature  of  the  plaama.  EquaUng  the  electrical  potenUal  onerg).  to  the 
IhormaJ  energy  wq  find  dial 


(8. 14) 


Solving  for  h  results  In 


ii 


(8. 15) 


The  quantity  h  defined  by  Eq.  (8. 15)  Is  called  the  Debye  ahleldlng  dietanee.  At  an  altitut 
o  a  few  hundred  kilometers  the  atmespherie  temperature  is  about  10®  -K  and  the  free 
e^tren  density  is  about  i0>®  eleetrens/m®.  If  we  assume  that  n^  in  the  above  equatio 
has  the  full  value  of  the  total  free  eleetron  density,  then  h  =  2  x  lo"®™. 


If  Uie  beam  has  a  length 
approximately 


L  and  a  radius  b,  the  average  distance  between  electrons  is 


»?  = 


(8.16) 


At  injection  bw  1  cm,  N  =  2  “i  v  j  »  ^ 

,  ,  ,  .  ’  ^  electrons,  and  L  =  1.2  x  10^  m.  (taking  the 

electron  s  velocity  equal  to  that  of  light,  and  a  pulse  duration  of  4  x  lo’^see) 


»?lnj  =  7  X  10"®m. 
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The  beam  will  blow  up,  boUj  in  cross  section  and  length,  until  tlie  average  distance 
between  electrons  is  roughly  the  Debye  distance  which  is  300  times  This  means 
that  b  will  Increase  to  about  3  m  and  L  to  about  4  x  lO^m. 

8. 3  BEAM  SPREAD  DUE  TO  MULTIPLE  SCATTERING  IN  THE  ATMOSPHERE 

Because  of  the  cumulative  effects  of  small  angle  scattering,  the  electron  beam  will 
spread  appreciably  during  its  passage  through^e  atmosphere.  The  importance  of 
this  spread  can  be  seen  from  the  values  of  as  a  function  of  altitude  h.  The 

angle  6  is  the  angle  between  the  electron  velocity  vector  and  the  vertical. 


h  {km) 


75 

0.34 

70 

0.58 

65 

0.85 

60 

»  1.0 

The  small-angle  approximations  used  in  computing  these  values  are  not  valid  for  0^  w  1, 
but  the  indication  is  clear  that  by  the  time  the  beam  gets  to  70  km  its  direction  is  no 
longer  well  defined. 


The  rapid  spread  of  particles  is  prevented  by  the  earth 's  magnetic  field  which  for  this 
calculation  will  be  assumed  vertical.  This  field  will  constrain  the  particles  so  that  they 
move  in  spirals,  thereby  limiting  their  horizontal  motion  between  collisions  to  the  diam¬ 
eter  of  the  spiral. 

An  approximate  calculation  will  be  carried  out  to  get  an  idea  of  how  important  the  effect 
is.  Numerical  calculations  indicate  that  during  most  of  the  passage  through  the  atmosphere, 
the  change  in  pitch  angle  per  cyclotron  revolution  is  smaller  than  the  instantaneous 
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pilch  angle.  Hence,  the  spread  of  particles  will  be  treated  as  a  diffusion  of  the  guiding 
centers.  At  each  altitude,  the  energy  will  be  assumed  monoenergetlc  and  equal  to  the 
average  energy  at  that  altitude.  The  gyration  radius  will  be  computed  at  each  altitude 

assuming  this  average  energy,  and  assuming  0  =\0  )  until  0  becomes  equal  to  60*, 
below  which  altitude  0  will  be  assumed  to  remain  at  60* .  The  time  a  particle  takes  In 
traversing  Ah  Is  equal  to  At  b  ojid  the  spread  In  the  beam  during  that  time 

win  be  equal  to 


Ar  =  At<Ar(0,n)>  =  At<Ar(h)>  (8.17) 

where  n  Is  the  number  of  atoms  cm  and  is  a  single-valued  function  of  h  ,  as  Is  0 
Thus,  the  problem  Involves  finding  an  expression  for  <Ar(h)>  which  Is  the  average 
change  per  unit  time  of  the  distance  between  the  guiding  center  of  a  particle  and  the 
centroid  of  the  electron  beam. 

With  the  geometry  shown  In  Fig.  8-2, 

r  =  distance  from  centroid  to  guiding  center  before  collision 
r'  =  distance  from  centroid  to  guiding  center  after  collision 
V'  =  projection  of  scattering  angle  on  horizontal  plane 
4  =  orbit  phase  angle 
a  =  cyclotron  radius 


mA  = 


2  2  2  2 
a  +  a  -  2a  cos  ^  =  2a  (1  -  coBip  ) 


J  ,2 

a  tjj 


r2. 


h^  -  2rh  cos  Jir  -  4  -  ~ 


2  2 

a  ^  -  2ra  ip  sin 


M) 


(8. 18) 

(8. 19) 

(8. 20) 
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Therefore, 


Fig.  8-3  Definition  of  Scatterli^  Angles 

gtn  _  sin  (b 
sin  ij  sin  O' 

sin  0  «  0  =  sin  n 

'  sin  6 ' 

cos  O'  =  cos  Tj  cos  6  +  sin  rj  sin  0  cos  0 


sin  0 '  -  (i  _  (cos  T/  cos  0  +  sin  rj  sin  0  cos 
sin  0 '  w  sin  0 


V 


sin  f 
sin  0 
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And  then, 


2 

1 

sln^  0 
2 


A-  =  _  sin  »  etn^  .  .  2  sin  6.  2. 

sin  e  - 2  (a  +  ra  cos  4) 


2r  sin  e 


(8.26) 


The  scattering  cross  section  for  electrons  in  air  will  be  taken  to  be 


-A-z  =  A  ^ 


4m^  o'*  Z?'*  sin'*  |  sin'*  | 

The  average  change  in  r  per  collision  is  given  by 

2ff 

[Arj  -  J  J  Ar((^,i;)  0(0  ,ij)  d0  dt;  sin  i; 

0  0 

Substituting  and  integrating  we  get 

|Arl=  4A  •  ^  (a^  +  ar  cosO  J 

sin  0  J  U 

**min  ^ 

IafI  =  — —  (a^  +  ar  cos  4)  /n  (— — \ 

'  '  r  sin^  0  V’^min/ 


(8.27) 


(8.28) 


(8.29) 


(8.  30) 
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Performing  the  average  over  orbital  phase  angle  {  , 


orbit 


'  r  0 


8»A 


\*^mln/ 


The  change  in  r  per  unit  time  is  then 


<Ar>  =  |Ar|  n  v 

^  _  Sir  Anv  .2 

<Ar>  =  - 2 —  ® 

r  sin  0 


(imli.) 


A(r  )  =  2r  Ar 


.  .  2.  „  dtr^l  _  16ir  Anv  _2  ._  /  2  \ 

Mr)  --ZzT^ 


sin^C 


(8.31) 


(8. 32) 


(8. 33) 


If  the  right  hand  side  is  expressed  as  a  function  of  t  by'  noting  that  h  =  h(t)  ,  the 
equation  can  then  be  integrated  numerically. 


dt 


=  f(t) 


2 


r 


(8.34) 


The  numerical  results  are  summarized  in  the  curve  of  Fig.  8-4. 
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0  10  20  30  40 

MEAN  SPREAD  (METERS) 

Fig.  8-4  The  Mean  Spread  Due  to  Multiple  Scattering  of  a  Narrow 

Beam  of  1  Mev  Electrons  Striking  the  Top  of  the  Atmosphere 
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8.4  RADAR  REFLECTION  FROM  AN  EXTENDED  REGION  OF  IONIZATION 

Combining  the  results  of  the  calculaUons  in  SecUons  8.2  and  8. 3,  the  region  of  ioni- 
saUon  is  found  to  have  the  foiiowing  dimensions.  In  alUtude,  it  wili  begin  about  70  km 
and  end  at  about  60  km;  in  iongitude,  it  wiii  have  a  width  of  about  100  m;  in  iaUtude,  it 
wiii  have  a  thickness  of  about  50  m.  As  computed  in  Section  8. 1,  there  wiii  be  about 
5  X  10  ion  pairs  distributed  throughout  this  voiumc.  If  wo  assume  that  the  distribuUon 
is  uniform,  then  there  will  be  an  Ion  pair  density  of  10^^  electrons/m®. 

We  wiii  assume  that  a  radar  staUon  can  dircci  its  beam  onto  this  "curtain”  of  ionizaUon, 
and  the  question  is  whether  or  not  there  will  be  a  detectable  echo.  (Fig.  8-5). 


Fig.  8-5  Relative  Position  of  Radar  Station  and  the  Region  of 
Ionization 


r 

r 
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The  lunount  of  ronccted  power  Is  determined  by  the  radar  transmitter's  frequency  and 
the  plasma  frequency  associated  with  the  lonlr.ed  region.  If  the  radar  frequency  Is  less 
than  the  plasma  frequency  there  Is  almost  total  renecUon,  and  If  It  Is  greater,  there  Is 
almost  no  reflecUon.  The  ma.ximum  frequency  of  the  radar  transmitter  Is  related  to  the 
Ion  pair  density  In  a  highly  reflective  Ionized  region  by 

2 

*'r  *  (8.35) 

where  Is  In  cps  and  N  In  Ion  pairs  per  m®.  In  our  case,  N=  10*^  Ion  paIrs/m® 
and  =  10  Mc/soc.  The  conclusion  Is  that  a  10  Mc/sec  radar  transmitter  would 
easily  detect  the  presence  of  the  electron  beam  In  the  upper  atmosphere. 

8.5  CONJUGATE-POINT  MEASUREMENT 

The  problem  of  radar  reflection  has  more  facets  than  Is  Indicated  by  the  simple  ealculatlon 
above.  These  addlUonal  complications  do  not  affect  the  basic  feasibility  of  the  method, 
but  must  be  examined  nonetheless. 

(1)  Due  to  the  Interference  phenomenon,  the  strength  of  the  reflected  signal 
may  depend  sensitively  on  the  angle  at  which  the  radar  beam  strikes  tlie 
ionized  region.  Normal  Incidence  will  produce  the  maximum  reflection. 

(2)  The  persistence  of  the  lonIzaUon  is  determined  by  the  rate  of  recombination. 

The  ionization  may  last  for  the  order  of  a  few  seconds  which  should  be  long 
enough  to  allow  the  radar  to  scan  and  lock  on. 

(3)  The  accuracy  with  which  the  position  of  the  ionized  region  can  be  located 
will  be  determined  by  the  angular  resolution  of  the  radar  transmitter  and  the 
distance  between  the  transmitter  and  the  reflector.  Because  of  the  llnc-of- 
sight  requirement,  a  particular  radar  station  will  have  a  maximum  range  of 
about  1,000  km  for  conjugate-point  detection.  If  the  beamwldth  of  the  10  Mc/scc 
radar  is  about  10° ,  then  at  the  maximum  range,  the  uncertainty  In  the 
conjugate-point  position  would  be  about  200  km.  This  poor  accuracy  is 
unacceptable.  The  accuracy  can  be  Improved  by  the  use  of  aircraft-borne 
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radar  so  that  the  transmitter  can  be  close  to  the  region  of  ionization.  If 
the  aircraft  is  within  10  km  of  the  position  of  the  ionization,  the  uncertainty 
in  the  measurement  would  be  about  2  km. 

(4)  Dyce  (Ref.  27 )  has  studied  the  problem  of  radar  reflection  from  auroral 
displays.  He  has  Indicated  that  it  would  be  possible  to  locate  the  poslUon  of 
the  ionized  region  with  ground-based  radar  with  greater  accuracy  by  (1) 
using  higher  frequency  radar  reflections  from  those  parts  of  the  ionized 
region  having  a  staUsUcally  higher  ion-i>alr  concentraUon,  and  (2)  making 
use  of  pairs  of  transmitters  and  receivers  as  a  radar  interferometer. 

(5)  Because  of  the  longitudinal  drift  of  the  electron  beam,  the  position  of  the 
ionized  region  will  be  sllghUy  east  of  the  magnetic  conjugate  to  the  injection 
point.  A  correction  can  be  made  to  the  data  to  take  account  of  this  fact.  In 

the  case  of  Injection  at  a  magnetic  latitude  of  45»,  the  eastward  displacement 
would  be  about  3  km. 

(6)  For  a  given  set  of  InjecUon  positions.  Table  8-1  gives  the  position  of  the 
conjugate  points  predicted  by  a  dipole  field  and  predicted  by  Welch's  512-term 
harmonic  expansion  of  the  geomagneUc  field.  It  is  the  purpose  of  this  com¬ 
parison  to  Indicate  the  magnitude  of  the  uncertainty  In  the  position  of  conjugate 
points  as  determined  by  two  different  geomagnetic  models. 
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Section  9 

SATELLITE  CHARGING 


0.1  INTRODUCTION 

The  mo.sl  riindamenlol  problem  which  must  be  examined  in  oitlor  to  demonstrate  the 
fcnsiblllty  of  artificial  Injection,  even  In  principle,  is  the  problem  of  satellite  charging. 

If  tile  satellite  were  moving  in  a  perfect  vacuum,  clcetrons  could  be  ejected  only  when 
the  potential  of  the  satellite  (due  to  the  residual  ix>sitivc  charge)  rose  to  a  value  equal  to 
the  electrons'  kinetic  energy.  At  this  point,  electrons  could  no  longer  leave  and  the 
useful  life  of  the  satellite  would  have  ended.  In  the  ease  of  pulsed  injcelion,  the  time 
for  a  satellite  to  rise  to  one  million  volts  would  be  I  psec,  for  d-e  injection,  it  would  be 
l/lO  sec.  Researchers  In  the  field  of  ionic  rocket  propiil.sion  have  studied  this  problem, 
and  Iiavc  found  the  only  suitable  solution  for  their  application  to  bo  the  simultaneous 
expulsion  of  positive  and  negative  charge  in  the  thrust  beam.  This  solution  is,  however, 
impractical  In  our  case  for  the  following  reason:  If  the  electrons  arc  not  to  be 
decelerated,  there  must  be  not  only  satellite  neiitriility,  .but,  also,  no  charge  separa¬ 
tion  in  the  ejected  ionized  bciain.  This  means  that  the  electrons  and  positive  ions  in  the 
beam  must  have  the  same  velocity.  If  the  positive  ions  are  protons  with  the  velocity 
of  a  onc-Mcv  electron,  then  the  proton  kinetic  cnergj'  is  2  Bev,  which  Is  obviously 
impractical.  Tlic  fact  is,  of  course,  the  satellite  is  not  moving  in  a  vacuum,  but  in  ,a 
rather  complicated  environment  made  up  of  the  earth's  magnetic  field  and  a  thermally 
agitated,  tenuous  atmosphere,  partially  neutral  and  partially  ionized.  Figures  9-1, 

9-2,  ,and  9-.3  show  the  orders  of  magnitude  associated  with  the  atmosphere  in  the  satel¬ 
lite  region.  The  ionized  component  of  the  atmosphere  is  a  source  of  free  eicetrona 
which  would  tend  to  neutralize  a  charging  satellite. 

The  effect  of  the  ambient  plasma  on  a  passive  satellite  has  been  studied  (Ref.  28),  and 
it  has  be;  n  shown  that  a  satellite  would  pick  up  more  electrons  than  positive  ions  (owing  to 
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Fig.  9-2  A  Normal,  Daytime,  Free-Electron  Distribution 
During  Sunspot  Maximum 
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the  difference  between  the  thermal  velocities  of  electrons  and  ions  at  Uie  same  tempera¬ 
ture)  and  would  rise  to  a  potential  between  10  and  100  volts.  For  a  rcasonable-sir.o 
satellite  ejecting  a  large  electron  current  and  charging  positively,  this  thermal  effect 
would  not  be  sufficient  to  maintain  neutrality.  However,  the  electric  field  of  the  charged 
satellite  would  draw  the  electrons  out  of  the  plasma  at  a  much  greater  rate.  The  ques¬ 
tion,  then,  is:  To  what  potential  will  a  satellite  rise  before  the  current  being  drawn  out 
of  the  ambient  plasma  is  sufficient  to  offset  the  ejeeted  current  and  limit  the  |)oicntlal? 

9.2  NEUTR.‘\LIZAT10N  BY  A  THERMAL  PLASMA 

It  would  be  possible  to  maintain  satellite  neutrality  by  tlic  electron  "tlicrmal  c.\ecss". 

The  size  of  tlic  satellite,  however,  would  have  to  be  very  large.  The  average  velocity 
of  particles  at  temperature  T  is 


v 


(9.1) 


where  k  is  Boltzmann's  constant,  and  T  is  in  degrees  Kelvin.  For  electrons  at 
1,000*  K  (0.1  cv) 


v  =  200  km /sec 

A  nitrogen  ion  at  the  same  temperature  has  a  thermal  velocity  of  1  km /.see.  Consequently, 
the  motion  of  the  plasma  ions  can  be  neglected.  Note  that  a  typical  satellite  velocity  is 
about  8  kin/sce;  ions  do  not  strike  the  satellite  but  arc  swept  up  by  it. 

The  thermal  electron  current  incident  on  the  surface  of  a  spherical  satellite  is 

I.j.h  =  47r  a^  n^  ev  (9.2) 

where  a  is  tlic  satellite's  radius  and  n^  is  the  plasma  electron  density.  For  the  two  oases 
of  pulsed  and  d-c  injection,  if  current  equilibrium  is  achieved,  the  following  satellite 
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radii  are  required.  (Tliroughout  the  study  of  satellite  charging,  it  will  be  assumed 
that  pulsed  injection  occurs  at  low  altitudes  (ng  ~10®  elcctrons/cm^  and  d-c  injee- 

n 

tion  at  high  altitudes  (n  ~  10  clcctrons/cm  3)  ,) 


*injectc<i "  *Th  "c  ” 

(amps)  (clcct/cm^)  (m) 

100  (pulsed)  10^  16 

10"^(d,c.)  10^  f) 


The  charging  problem  can  therefore  bo  solvetl  in  the  last  resort  by  using  an  Echo-typo 
satellite  shell. 


In  the  case  of  pulsed  injection  with  a  pulse  repetition  rate  of  about  one  per  second,  a 
normal-size  satellite  charged  as  high  as  a  million  volts  would  Ixj  neutralized  completely 
by  the  thermal  electrons  in  the  time  between  pulses. 

9.3  THE  EFFECT  OF  A  PLASMA  ON  A  CHARGED  SATELLITE 

9.3. 1  Description  of  the  Model  to  be  Studicil 


Let  us  take  a  simple  model  in  which  the  the  following  assumptions  are  made: 

•  The  satellite  is  a  spherical  conductor  with  a  radius  a  ,  and  is  ejecting  a 
constant  electron  current  I^. 

•  The  effect  of  the  earth's  magnetic  field  is  neglected,  as  well  as  any  magnetic 
effects  due  to  the  ejected  current  or  the  flow  of  plasma  electrons. 

•  Initially,  the  plasma  is  neutral  and  has  a  constant  electron  number  density  n 

•  The  satellite  and  the  ambient  ions  are  at  rest,  and  the  thermal  motion  of  the 
plasma  electrons  is  neglected. 

•  In  describing  the  effect  of  an  electric  field  on  the  plasma  electrons,  it  will  be 
assumed  that  it  is  possible  to  assign  a  transient  conductivity  cr^  to  the 
medium;  that  is,  a  conductivity  which  is  a  function  of  time  but  which  approaches 
asymptotically  a  constant  value  cr^  characteristic  of  a  steady-state  current  flow 
in  the  plasma. 
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Splizcr  (Rof.  12)  has  descrilKicI  the  problems  Involved  In  defining  a  plasma 
conduelivlly  and  gives  the  steady-stale  values  for  various  models. 

9.3.2  Transient  Conductivity 


The  equation  of  motion  of  an  electron  in  the  plasma  is 


m 


dt 


c  E 


or  V 


(9.3) 


uhcic  the  effect  of  collisions  within  the  plasma  has  been  represented  bj'  a  retarding 
force  proportional  to  the  electron's  velocity.  The  constant  a  can  bo  evaluated  by 
noting  that  in  the  steady  state 


!§  =  or  V 


The  current  density  is  given  by 


j  =  e  n  V 
c 


Combining  Eqs.  (9.4)  and  (9.5)  we  obtain 


^  P  -► 


(9.4) 


(9.5) 


(9.0) 


Therefore, 


2 

n  e 
e 


O’ 

o 


(9.7) 
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If  the  electric  flel.l  hoes  not  ory  a|)|)recl:,l.ly  over  the  tllslancc  between  collision  ccnlcrs 
or  In  times  com,.:,,:, We  will,  the  eolllslon  period,  then  Eq.  10.3)  mny  be  lnleBr.nled 


_ o_ 

n  c 
c 


1  -  exp 


-&<)■ 


E 


(9.8) 


where  it  hv.s  been  assuinecl  llmi  the  Inllliil  velocity  is  zero.  Applying  Eq.  (9.r,) 


j  ^  0^ 


1  -  exp  - 1 


n  e 
c 

m  O’ 


2  \  1 


E 


(9.9) 


The  transient  conductivity  is  then 


(T=  a 


1  -  exp  - 


n  o 
e 

m  a 


.2v 


(9.10) 


For  the  simple  plasma  wc  have  assumed,  is  given  by  (Ref.  12) 


O’  = 


y  rj,. 

E 


,3/2 


o  (:}8)  Z  in  A 


mho/m 


(9.11) 


where  T  Is  the  pl.tsma  temperature.  Z  is  the  Ionic  charge,  and  the  values  of  In  A 
and  are  given  in  Table  0-1.  The  eonduetivity  is  shown  in  Fig,  0-4,  assuming 

the  phase, a  is  singly  Ionised  (Z  .  1)  .  The  remarkable  feature  Is  the  eonstaney  of  the 
conductivity  over  the  large  range  of  altitudes. 


The  v.ilues  of  the  quantity  (n^eVmo^)  appearing  in  the  exponent  In  Eq,  (9. 10)  a 

given  in  the  tabic  below. 
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Injection 
pulsed 
d.  c. 


"e(°‘“Vcm») 


X{m. ) 


In  the  pulsed  injection,  the  pulse  duration  is  4  x  lO"®  sec;  the  exponential  in 
Eq.  (9. 10)  can  therefore  be  replaced  by 


giving 


(9. 12) 


(9. 13) 


This  equation  shows  that  the  plasma  current  is  limited  by  the  electron's  inertia 
rather  than  by  collisions.  In  the  d-c  case,  the  value  of  the  exponent  Indicates  that 
the  time  for  the  steady  state  to  be  achieved  is  about  10  sec. 


The  quantity 
values  of  the 


f4)" 


1  equal  to  the  plasma  collision  frequency  which  leads  to  the 
mean  free  path  for  electrons  X  appearing  in  the  table  above. 


9. 3. 3  The  Satellite  Potential 


We  are  now  in  a  position  to  write  an  equation  for  the  potential  on  the  satellite  as  a 
function  of  time. 


V 

B 


(9.14) 
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The  capacity  of  the  satellite  Is 


C  =  4iff  a 
K  o 


(9. 15) 


(There  lx  some  .irbllrnrlncss  In  the  permllUvIly  sssigned  to  n  pUismn.  but  wllh  the 
eondMlIrli).  we  have  used,  the  permittivity  will  have  Us  vacuum  value  t  ). 

O' 

The  satellite  charge  is 


2 

*>8  j  J  (a.  t')df 


(9. 16) 


From  the  definition  of  the  transient  conductivity 

J(a,  t)  =  e^Efa,  t) 


=  a 


Is. 


t  4  jrc  a2 
o 


(9. 17) 


Therefore , 


%  = 


^  I 

1“  / 

V  /V 


(9. 18) 


Differentiating  both  sides  we  have 


do  (7 

— §  a.  t 

dt  %  ~ 


(9.19) 


The  solution  of  this  equation  which  satisfies  the  initial  condition  q^(0)  =  0  is 
[exp  -  J  dl’j  j‘dl.  exp  dl" 

®  -1  O  n 


Qg  =  |exp 


(9.  20) 
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For  tho  case  of  pulsed  Iiilectlon  <7^  Is  given  by  Eq.  (9. 13). 

"s  = -  (i^)  4  ]  /  <»'  -p  4 


(9.  21) 


The  quantity 


(^)» 


the  square  of  the  plasma  frequency  .  For  an  electron 


density  of  10  electrons  per  cm^ ,  w  ^  =  5. 65  x  10^  rad/sec.  Equation  (9. 21) 

can  be  Integrated  approximately  to  give 

2 


"'s  ' 


(0.3)  oxp^-  -f  ^1  -  oxp  0,60.2 


(9.22) 


In  two  limiting  cases,  q  becomes 

8 

"s  ■  V‘(*  ■ 


I  I 

P 


2  2 

a,;t"  <«  1. 


2.2  »  1. 
w  t 
P 


(9.  23) 


From  Eqs.  (9. 14)  and  (9. 22),  the  satellite  potential  can  be  computed  and  Is  shown 
In  Fig.  9-5.  The  potential  Initially  rises  as  If  there  vfrcre  no  plasma  present.  It 
Is  during  this  time  that  the  electron's  Inertia  Is  being  overcome.  Then,  after  about 
10  seconds  have  elapsed,  a  significant  number  of  electrons  start  pouring  In  and 
the  potential  flattens  off  and  then  falls.  At  the  end  of  the  pulse, q  decays  at  a 
much  greater  rate 


(j) 


=  qg(T)  exp 


(t2  -  t2 


) 


t  >  T 


(9.  24) 
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During  the  Electron  Injection 


where  T  Is  the  pulse  (hsr.itlon,  Tl»e  pulsed  accelerator  described  In  Seetinn  .1.3 
has  an  Intrinsic  enecKj'  spread  of  about  *1  percent  of  the  Injection  energj-  or  1  x  10 

volts.  Tlds  means  that  the  effect  of  the  satellite  cluarglnR  would  bo  osson- 
tlnlly  ncRliRlble. 


-1 


In  Ujo  ease  of  d-c  Injection, 


“'j,  =  5.65  X  10-’  rad/scc’'  („  =  10^  clcetrons/cm‘* 


UshiR  Eq.  (n.  22)  to  compute  the  tninslent  charRlnR  of  the  .satellite  at  the  time  of  the 
Initiation  of  Injection.  It  Is  found  that  the  satellite  potential  rises  to  a  maximum  of 
10  volts.  After  alxnit  10  .seconds  Im  e  eh, p.ses.  ean  be  replaeed  by  o  .  „„d  the 

imtcatlal  achieves  Its  stc.ady-.statc  value  of  a  fraotlou  of  a  volt. 

9. 3.  -I  Validity  of  the  Motlcl 


For  U,o  model  we  have  taken.  It  Is  clear  that  satellite  charging  Is  not  a  serious 
problem.  However,  the  question  of  the  validity  of  the  assumptions  Involved  In  making 
the  calculations  must  be  examined.  It  will  not  be  possible  to  rigorously  justify  the 
model,  but  a  number  of  arguments  can  be  made  which  Indicate  that  the  m.ixlmum 
satellite  potential  will  not  be  vastly  different  from  that  computed. 

Geomagnetic  field  makes  conductivity  anlsotronlc.  The  presence  of  the  earth's 
magnetic  field  has  the  effect  of  m.Tklng  the  conductivity  no  longer  isotropic.  In  the 
direction  parallel  to  the  magnetic  field,  the  conductivity  will  be  equal  to  the  value 
we  have  used.  Perpendicular  to  the  magnetic  field,  however,  the  conductivity  will 
be  many  orders  of  magnitude  smaller.  Our  model  can  be  applied  to  this  situation 
if  we  allow  current  to  flow  onto  the  satellite  only  along  the  magnetic  field  lines, 

setting  the  conductivity  in  other  directions  equal  to  /.ero.  In  this  case,  Eq.  (9. 16) 
must  be  replaced  by 

TT 

2  ‘ 


IT 

rT 


2  r2  2  r 

Tra  J  do  cos  0  sin  0  j  dt'  j  (a,  t'). 

fl 


(9.  25) 


o 


9-15 


LOCKHEED  AIRCRAFT  CORPORATION 


MISSILES  ana  SPACE  DIVISION 


2-44-61-1 


Then, 


q 


s 


-  4/3)ra^ 


(9. 26) 


whore  j(a.t)  remains  the  current  density  at  the  satellite  surface  with  no  magnetle  field 
present  The  solution  of  Eq.  (9.26)  is  the  same  as  Eq.  (9. 22)  with  w  ^  replaced  by 
1/3  Wp  .  This  will  Increase  the  maximum  satellite  potential  by  a  factor  of  3^^^ 

Pjasma  density  remains  nearly  constanc.  In  relating  T.  v".  and  E.  it  was  assumed 
that  the  plasma  electron  density  remains  constant.  To  make  this  assumption  rea¬ 
sonable.  It  can  be  shown  that  If  at  some  time  the  plasma  eharge  density  p  in  Uie 
vicinity  of  the  satellite  differs  from  zero,  the  plasma  will  return  to  neutrality  in  a 
time  which  is  short  compared  with  the  time  for  any  significant  change  of  Uie  satellite 
potential  to  occur.  From  Maxwell’s  equations  and  the  condition  of  charge 
conservation: 


(9. 27) 


^  at 

J  and  E  are  related  by 


Solving  for  p  from  these  three  equations 


and  Integrating  we  get 


(9.28) 


(9.29) 


(9.30) 


(9.31) 
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Ifwu  take  tlio  value  of  from  Kq.  (9.13),  then 

2 

tii 

P  =  PqCxp--^  t  results  (9,32) 

The  time  Interval  for  p  to  fall  to  e''  of  Its  Initial  value  Is 


2  -8 

^  -1  ~  ~  =  2.5  X  10  sec  (9.33) 

e  p 

If  we  look  at  Fig.  9-5,  we  see  that  the  time  for  the  satellite  potential  to  fall 
to  e"  of  Its  mioclmum  value  Is  about  lO’"^  sec,  or  four  times  Uie  plasma  neu¬ 
tralization  time  constant.  This  argument  does  not  hold  In  the  time  Interval  In 
which  the  satellite  Is  rising  to  Its  ma.xlmum.  However,  In  this  Interval,  the  rate 
of  satellite  charging  Is  essentially  ^hat  with  no  plasma  present,  and  the  assumption 

of  a  constant  n^  (and  a  constant  E  ,  as  discussed  In  the  next  paragraph)  Is  not 
required. 

The  effect  of  a  varying  electric  field.  In  order  to  be  able  to  Integrate  the  equation 

of  motion  of  a  plasma  electron  (Eq.  (9.3)),  the  electric  field  had  to  be  constant 

to  the  degree  that  the  energy  gained  between  collisions  was  nearly  constant.  From 

the  computed  values  of  the  mean  free  paths  and  collision  frequencies.  It  is  clear 

that  this  assumption  is  invalid.  Let  us  calculate  the  satellite  charge  again  without 

this  restriction  on  the  eleetrie  field,  but  taking  n  constant  and  neglecting  collisions 
altogether. 


(9.  34) 


(9.35) 


.1  = 


e  n  V 
n 


(9.36) 
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%  (1) 

E(a,  t)  =  - -  (9,37) 

4  »c  a 
o 


Differentiating  Eq.  (9.35)  twice  and  combining  the  above  equations  we  find  that 


.2  2 

^  ®  "e  2 

■“2“  "  -irr  Is "  -“p",  («•“) 


The  solution  of  this  equation  Is 


q 


s 


sin  u)  t 
P 


(9.39) 


With  the  modified  assumptions,  then,  q^  oscillates.  The  important  thing,  however, 

is  that  the  mxxlmum  value  of  q  leads  to  a  maximum  satellite  potential  of 
4  ® 

2x10  volts  for  pulsed  Injection,  and  20  volts  for  d-c  injection.  These  values  arc 
approximately  equal  to  those  computed  before. 


In  conclusion,  the  results  of  this  analysis  indicate  that  the  maximum  potential  to 
which  the  satellite  will  rise  is  about  3  x  10**  volts  in  the  pulsed  case,  and  about 
30  volts  in  the  d-c  case.  Therefore,  satellite  charging  will  not  seriously  affect  the 
monoenergetic  property  of  the  electron  beam  for  the  accelerators  considered  In 
this  study. 
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Section  10 

DIAMAGNETIC  EFFECT  OF  THE  INJECTED  ELECTRONS 
10.1  INTRODUCTION 


It  will  be  necessary  to  compute  the  magneUc  field  associated  with  Uie  injected  electron 
current  to  show  that  it  is  small  compared  with  the  earth's  field  In  the  region  In  which  the 
electrons  are  moving.  If  this  were  not  the  case,  the  electrons  would  see  a  modified  field, 
and  It  would  be  much  more  difficult  to  correlate  the  results  of  the  various  experiments 
with  similar  phenomena  occurring  in  the  unperturbed  field.  In  the  worst  case,  the  earUi's 
field  could  be  reduced  (the  electron  current  produces  a  diamagnetic  field)  to  the  point 
where  the  electrons  would  no  longer  be  magneUcally  contained. 

10.2  PULSED  INJECTION 

The  cyclotron  period  of  a  high-energy  electron  In  a  magnetic  field  is 


T  = 

c  c  Be 


(10. 1) 


where  p  Is  the  electron's  momentum.  For  a  1-Mev  electron  circulating  near  the 
surface  of  the  earth,  the  cyclotron  period  is  nearly  equal  to  the  duration  of  the  Injection 
pulse  (4  X  10  sec)  ,  and  becomes  much  longer  far  from  the  earth.  Therefore,  the 
pulse  of  electrons  can  be  represented  by  a  single  current  loop  carrying  100  amps. 

From  Smythe  (Ref.  1),  the  components  of  the  magnetic  field  of  a  current  loop  are 
given  by  the  following  equations,  where  the  coordinate  system  is  defined  In  Fig,  10-1. 
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Fig.  10-1  Coordinate  System  of  a  Current  Loop 


B 

^  V . 

1 

r.  .  .  ,2  1 

z 

2jr  1 

(a  +  p)^  + 

172 

- 2 - T 

(a  -  pr  +  z  J 

(10.2) 

B 

P  I 

_  o 

z 

\  .  a^  +  p2  2 

P 

27r 

P 

[(a  +  p)2  + 

ji72 

-K  +  - £- — - — £_  p 

(a  -  p)2  +  z^ 

(10.3) 

K 

rV2 

do 

0 

(l  -  sln^  e 

^2 

(10.4) 

/■V2  /  1/2 

E  d0(l-k2  8in2 


10-2 


LOCKHEED  AIRCRAFT  CORPORATION 


MISSILES  and  SPACE  DIVISION 


2-44-61-1 


4ap 


(a  +  p)^  +  7?" 


(10.6) 


To  compute  the  field  In  the  vicinity  of  the  current  filament  itself,  wo  sot 


p  =  a  -  c 
f/a  «1 
z  «a,  p 


Then, 


(10.7) 


Z  «  € 

«  1  - 

1  /  2  ,  2. 

— 2  +  z  ) 

4a 

(10.8) 

K  »  Xn 

8a 

(10.9) 

E  1 

(10. 10) 

Keeping  only  the  dominant  terms,  we  find  that  B  becomes 

z 


R  =  fi-  +  -L  « 

®z  27r  [^2  2a  r  J 


,1/2 


(10.11) 


/  2  2\  ' 

where  r  =  +  z  j  ,  which  is  the  radial  distance  from  the  center  of  the  current 

filament  (see  Fig.  10-2). 
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z 


Fig.  10-2  Cross  Section  of  a  Current  Loop 

•nie  first  terra  In  Eq.  (lo.  11)  is  siraply  the  ■'straight  wire"  raagneUe  field,  and  the 

seeond  term  Is  the  effect  of  the  curvature  of  the  ■Vlre."  We  will  neglect  the  first  terra 

beeause  It  produces  a  raagnetle  "pineh."  but  dees  not  result  In  a  diamagnetic  field;  the 
second  term  does. 


fi 


dia 


8a 

47ra  r 


(10.12) 


As  two,  diverges.  However ,  the  current  filament  must  not  be  taken  as 

infinitely  thin,  but  as  having  some  finite  spaUal  distribution.  We  will  define  a  current 

ensity  ,v,  ,,,  ^ 

radius  of  the  current  filament. 


j  = 


(10. 13) 
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The  diamagnetic  field  is  then 


=  llm  -  ~  -1- 


'dla 


6-,0  *2 


.r-6 


r*  “  /•"  /2.2 

j  pdpj  (l(pl.n  - 2 — —  2r  P  cos  ij>) 


61a 


+  /  p  rip  /  V  f  „(r^  +  -  2r  P  cos 


r+6 


64a 


(10. 14) 


Integrating  and  keeping  only  the  dominant  term  we  get 


B 


dia 


4jra  ^  b 


(10.15) 


The  diamagnetic  field  ,  then,  Is  essentially  constant  In  the  vicinity  of  the  current 
filament.  The  radius  of  curvalnre  a  Is  related  to  Ihe  magnetic  field  by 


Therefore , 


(10. 16) 


B u  I 

_  lo_  ce  8a 

B  4;r  cp  "  b  (10.17) 


In  the  case  of  a  1-Mev  electron  clrculaUng  near  the  surface  of  fte  earth,  the  para 
meters  In  Eq.  (10-17)  have  the  following  values  at  a  magnetic  latitude  of  SO". 

I  =  100  amps 

B  =  0.34  gauss  =  3,4  x  10“®  webers/m^ 
a  =  140  m 

cp  =  1.4  Mev  =  2.24  X  10~^^  joules 
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And  if  b  =  1  cm,  then 


B 


dia  ,  -  -2 

Q—  =  1.6  X  10 


'n»e  effect  of  the  diamagnetic  field  Is.  therefore,  negligible. 

>0-3  DC  INJECTION 

In  the  cane  of  <l-e  Injection.  U.e  argument  of  the  last  section  does  not  applj.  beeanse 
the  electron  current  la  distributed  over  .  large  volume  of  space  and  eannet  be  ropre- 

sentod  by  some  simple  configuration  such  as  a  loop  or  a  solenoid.  We  . . . 

assume  that  the  densHy  of  Injected  electrons  is  high  enough  so  that  we  can  assign  a ' 

volume  magnetisation  to  the  region  occupied  by  the  electrons.  Associated  with  each 
electron  is  a  magnetic  moment 


M  =  jre  n  a 
c 


2  2 
amps-m 


(10. 18) 


where  Is  the  cyclotron  frequency  and  a  Is  the  cyclotron  radius.  Taking  the 
olcotron’s  velocity  to  be  equal  to  c  .  we  find  that 


I'  = 


<v 

c 
2  n 


c 

27ra 


a  = 


Be 


(10.19) 


Therefore, 


M  =  §1 

2B 


(10.20) 
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I 


t 

I. 


i . 


f 


The  electron  number  density  is  n  ,  and  tiie  magnetic  moment  per  unit  volume  is 

0 

“  "e  2B  amp/m  (10.21) 

The  quantity  M  Is  called  the  magnetization,  and  Is  the  source  of  a  diamagnetic  field 


“dto  °  "o  "e  2B 


(10.22) 


We  must  now  digress  In  order  to  compute  the  value  of  n^.  In  the  experiment  dcs- 
crlbcxl  In  Section  6.1,  It  is  proposed  to  inject  electrons  Into  a  volume,  the  cross  section 
of  which  Is  shown  in  Fig.  G-1. 

The  shaded  region  In  this  figure  Is  defined  approximately  by  the  two  field  lines  which 
cross  the  magnetic  equator  at  (r^  -  a^)  and  (r^  +  a^)  ,  where  a^  is 


\  ^  cB(r  =  r*^  ,  tl)=0)  efi"  (10.23) 

e  e 

For  a  dipole  field,  the  equation  of  a  field  line  is  gdven  by 

2 

r  =  r  cos  Ip  (10.24) 

0 

The  volume  of  the  "magnnetic  shell"  is  represented  by 


i 
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V  w  3 


•  «5i  a  r 
e  o 


a  «r 
e  e 


(10. 26) 


TTie  electron  number  density  Is  approximately 


I  T 

n  -  X  _D 

e  e  V 


(10. 27) 


D  Is  11.0  Umc  for  the  electrons  to  drift  longltodlnally  once  around  the  earth  which  Is 
opproxlmately  the  length  of  the  InJeeUon  Ume.  From  Eq.  3.  Id. 


P 

m  ^  4  TT  © 

D  3  a  c 
e 


(10. 28) 


Going  back  to  the  expression  for  the  diamagnetic  field  [Eq.  (10.22)]  ,  and  suhsUtutlng 
from  Eqs.  (10.26),  (10.27),  and  (10.28) 


dia  cp  p  ( 

-r-eo"e^  =  (0.2,Poi,fy 


(10.29) 


This  ratio  has  its  maximum  value  at  the  equator  (B  =  B  ). 


B 


dia  _  .  e 

B  -  F 


For  one  Mev  electrons  and  =  io~^ 


amps 


(10.30) 


B 


dia 


B 


=  5  X  10'® 


The  diamagnetic  effect  is  completely  negligible. 
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Using  (ho  above  fonnnlas.  It  might  bo  Intorosling  If  wo  oom,x,tc  „„ 
of  Iho  nohlovnble  oloetron  density  In  a  filled  magnolle  shell. 


appro.xinwtc  value 


n 


0 


(0.4) 


For  a  dipole  field, 

B  =  (9. 1  X  10^®)  r  ^  Webers/  2 
®  m 

Taking  the  case  where  fg  =  2  earth  radii, 


(10.31) 


n  =7  electrons/  3 
m 

which  gives  a  Dux  of  2  X  lo®  electrons  per  cm^’  per  sec. 
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Section  11 

CONCLUSIONS  AND  RECOMMENDATIONS 


11.1  CONCLUSIONS 


Three  types  of  experiments  in  which  electrons  are  injected  artificially  into  the  earth's 
magnetic  field  by  satellite-borne  accelerators  have  been  considered  and  appear  to  be 
feasible.  Each  of  these  cxpcilments  has  considerable  scientific  merit  and  militarj- 
usefulness.  The  experiments  arc: 

(1)  Integral  Invariant  surfaces  can  be  covered  by  fast  electrons  which  have  been 
injected  by  a  d-e  accelerator  operating  during  the  passage  of  the  satellite 
through  the  surface.  The  longitudinal  drift  of  the  electrons  serves  to  close 
the  shell  and  so  can  provide  a  measure  of  the  longitudinal  drift  time.  In 
addition  to  this,  tlie  electron  life-times,  the  mechanisms  determining  the 
life-times,  and  diurnal  and  magnetic  storm  effects  on  the  shape  of  the  surface 
can  be  observed.  It  appears  that  inclined  elliptical  orbits  with  electrons 
injected  near  apogee  are  probably  most  suitable  for  both  injection  and  sub¬ 
sequent  observation  of  the  electrons. 

(2)  The  bounce  time  of  the  electrons  Injected  by  a  pulsed  accelerator  can  be 
measured  by  detection  of  the  synehroton  radiation  with  a  radio  receiver 
carried  in  the  injecting  satellite.  Injection  would  be  made  at  magnetic  pitch 
angles  of  90  deg  with  a  pulse  duration  of  four  nsec.  After  one  or  more 
bounces,  the  electrons  will  be  mirrored  in  the  vicinity  of  the  satellite  and 
the  time  between  injection  and  subsequent  mirroring  can  be  measured.  This 
experiment  will  measure  the  average  value  of  the  magnetic  field  along  the 
field  line  on  which  injection  is  made  and  will  permit  comparison  of  this 
average  with  that  given  by  various  geomagnetic  field  models.  The  variation 
of  the  magnetic  field  produced  by  diurnal  effects  and  by  magnetic  storms  can 
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bo  observed  as  can  the  imucimum  laUtude  at  which  closed  lines  occur.  A 
polar  orbit  satellite  should  be  used  for  this  experiment. 

(3)  GeomagneUc  conjugate  points  can  be  measured  by  Injecting  electrons  along 
a  magnetic  field  line  and  observing  the  points  at  which  the  electrons  penetrate 
the  atmosphere  In  opposite  magnetic  hemispheres.  The  electrons  would  bo 
Injected  from  a  pulsed  accelerator,  carried  by  a  polar  orbit  satellite.  The 

establishment  of  conjugate  points  Is  particularly  Important  to  mllitaiy  appli¬ 
cations. 


11.2  RECOMMENDATIONS 

since  Iheec  experimenUi  appear  to  be  feaelble  and  deelroble,  U  Is  recommended  that 
further  effort  be  put  Into  opUmUIng  Ihc  physio's  design  of  the  experiments  and  Into  the 
development  of  an  engineering  program  plan.  The  physic's  design  of  the  experiments 
Should  receive  further  attention  In  the  following  areas; 

(1)  In  order  to  optimize  the  Injeetion  and  subsequent  deteetion  of  eleetrons  by  a 
satellite  eariylng  a  d-e  aeeelerator,  the  motion  of  the  satellite  In  the  magnetie 
field  of  the  earth  should  be  eonsldered  In  detail.  This  eonslderation  should 
include  the  Inclination,  eeeentrielftr,  precession,  and  decay  of  the  satellite  orbit 

and  the  rotation  of  Ihe  earth  and  Its  mngnetle  field.  The  study  will  require 
maehlne  ealeulation. 

(2)  Measurements  of  the  natural  eleetron  fiux  whleh  will  eontribute  a  baekground 
for  eleetron  deteetors  should  be  pursued.  These  measurements  are  eurrently 
in  progress  at  a  number  of  ageneles,  but  some  additional  measurement  are 
required.  In  partieular.  the  flux  of  eleetrons  having  energies  near  one  mev 
should  be  established  at  different  earth  radii. 

(3)  The  teehnique  of  deteeting  the  synchroton  radiation  from  a  vehiele-borne  radio 
receiver  should  be  explored  further.  It  may  be  possible  to  observe  the  natural 
synchroton  radiation  from  Van  Allen  eleetrons  or  from  auroral  eleetron  fluxes 
with  receivers  carried  in  sounding  rockets  or  satellites. 
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(4)  The  detection  of  atmospheric  ionization  in  Urn  conjugate  point  experiment 
should  be  further  investigated.  Fruitful  areas  for  investigation  include 
looking  nt  auroras  witii  radar  and  measuring  the  electron  flux  which  produces 
visible  auroras.  It  may  be  possible  to  do  the  conjugate  point  experiment  by 
optical  means  rather  than  by  radar. 

(5)  An  engineering  program  plan  and  cost  estimate  should  be  developed  at  Uic 
completion  of  the  additional  physic's  study.  The  program  plan  should  bo 
developed  after  careful  consideration  of  special  equipment  which  will  be 
required  both  in  the  satellite  and  on  the  ground.  Development  of  special 
equipment  such  as  nccclemtors,  low-noise  radio  receivers, and  vehicle 
memoiy  systems  will  probably  be  required. 
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No.  of  Cys. 

HEADQUARTERS  USAF 
1  Hq  USAF  (AFDAP).  Wash  25.  DC 

1  Hq  USAF  (AFORQ).  Wash  25.  DC 

1  Hq  USAF  (AFDRT-GW/1).  Wash  25.  DC 
1  Hq  USAF  lAFClN).  Wash  25.  DC 
1  '  Hq  USAF  (AFTAC).  Wash  25.  DC 

ATIC.  Wright-Pattorson  AFB.  Ohio 
1  (AFCIN  4F2A.  Maj  Pearce) 

1  (AFCIN  4DIA.  Library) 

I  AFOAR  (RRN.  Maj  Mtinyon).  Bldg  T-D.  Wash  25.  DC 


MAJOR  AIR  COMMANDS 
AFSC.  Andrews  AFB.  Wash  25.  DC 
1  (SCR) 

1  (SCR-2.  Capt  Ray  Berrier) 

1  Office  of  Aerospace.  Research.  AFOSR.  (SREC).  Wash  25.  DC 

AFSC  ORGANIZATIONS 

Office  of  Aerospace  Research.  Hq  AFCRL.  L.  G.  Hanscom  Fid. 
Bedford.  Mass 

1  (Tech  Library  for;  Mr.  Horowitz) 

1  (CRZCT,  H.  P.  Gavin) 

1  BSD  (Tech  Library).  AF  Unit  Post  Office.  Los  Angeles  45.  Calif 

1  SSD  (TDC  61-372-2).  AF  Unit  Post  Office.  Los  Angeles  45.  Calif 

1  ASD  (WWAD).  Wright-Patterson  AFB.  Ohio 

1  ESD  (Tech  Library).  Hanscom  Fid.  Bedford.  Mass 

1  AFMTC  (MTASI).  Patrick  AFB.  Fla 

1  Aeronautical  Research  Lab  (RRLO.  Mr.  Cady). 

Wright-Patterson  AFB.  Ohio 

1  AEDC  (AEOI).  Arnold  AFS,  Tcnn 

I  APGC  (PGTRI).  Eglin  AFB.  Fla 

1  RADC  (RCOII,-?.).  Grifflss  AFB.  NY 

AFSWe.  Kirtland  AFB,  NMcx 

1  (SV/NH) 

60  (SWOI) 

1  (SWRPA) 

1  (SWRPS,  Maj  Nadler) 
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1  {S\7RPT.  dipt  Welch) 

I  (SWRIN..  Capt  O'Brien) 

1  (SWRJ) 

OTHER  AIR  FORCE  AGENCIES 

1  Diicctor.  USAF  Project  RAND,  via:  Air  Force  Liaison  Office. 
The  RAND  Corporation  (RAND  L<ibrary  for:  Dr.  A,  Latter), 
1700  Main  Street,  Santa  Mimica,  Calif 
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1 
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ARMY  ACTIVITIES 

Chief  of  Research  and  Development,  Department  of  the  Army, 
{Special  Weapons  and  Air  Defense  Division),  Wash  25,  DC 

ARGMA  (  Tech  L.ibrary).  Huntsville,  Ala 

Director,  Ballistic  Research  Laboratories  (Library  for:  Mr,  E  O 
Baicy),  Aberdeen  Proving  Ground,  Md 

Director,  Evans  Signal  Laboratory  (Weapons  Effects  Section) 
Belmar,  NJ 

Commanding  Officer,  US  Army  Signal  Research  and  Development 
Laboratory  (Tech  Library  for;  Dr.  W.  McAfee),  Ft  Monmouth,  NJ 

NAVY  ACTIVITIES 

Chief  of  Naval  Operations,  Department  of  the  Navy  {OP-36) 

Wash  25,  DC 

Commanding  Officer,  Naval  Research  Laboratory  (Tech  Librarv) 
Wash  25,  DC  ^  i-mraryj. 

Commanding  Officer,  Naval  Radiological  Defense  Laboratory, 
(Technical  Information  Division),  San  Francisco  24,  Calif 

Office  of  Naval  Research  (Dr.  W.  J.  Thaler),  Wash  25,  DC 


OTHER  DOD  ACTIVITIES 

Director,  Weapon  Systems  Evaluation  Group,  Room  2E1006, 

The  Pentagon,  (Dr.  Gerald  Rosen),  V/ash  25,  DC 

Chief,  Defense  Atomic  Support  Agency  (Document  Library  Branch 
for:  Maj  Ralph  Pennington),  Wash  25,  DC 

Commander,  Field  Command,  Defense  Atomic  Support  Aeenev 
(FCAG3),  Sandia  Base,  NMex 

ASTIA  (TIPDR),  Arlington  Hall  Sta,  Arlington  12,  Va 

Director,  Advanced  Research  Projects  Agency,  Department  of 
Defense,  The  Pentagon  (Lt  Col  Roy  Weidler),  Wash  25,  DC 


AEC  ACTIVITIES 

US  Atomic  Energy  Commission  (Technical  Reports  Library 
Mrs.  J.  O'Leary  for  DMA),  Wash  25,  DC 
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I  President,  Sandia  Corporation  (Document  Control  Divisior  for: 

Dr.  Torn  Cook),  Sandia  Base,  NMex 

I  Chief,  Technical  Information  Service  Extension,  US  Atomic  Energy- 

Commission,  Box  62,  Oak  Ridge,  Tenn 

I  Director,  University  of  California  Lawrence  Radiation  I.nboratory, 

(Technical  Information  Division,  Mr.  Clovis  Craig  for;  Mr.  N.  C. 
Christofilos),  P.  O.  Box  808,  Livermore,  Calif 

I  Director,  Los  Alamos  Scientific  Laboratory  (Helen  Redman,  Report 

Library),  P.  O,  Box  1663,  Los  Alamos,  NMex 

OTHER 

1  Institute  for  Defense  Analysis,  Room  2B257,  The  Pentagon, 

(Dr.  \V.  Bradley),  Wash  25,  DC 

1  General  Electric  TEMPO  (Dr,  R.  W.  Hendrick),  735  State  Street, 
Santa  Barbara,  Calif  Contract  AF  29(601)-4194 

1  Lockheed  Missile  &  Space  Division  (Dr.  Roland  E.  Meyerotl), 

3251  Hanover  St.  ,  Palo  Alto,  Calif  Contract  AF  29(601)-4 1 36 

i  E.  H.  Plesset  Associates  (Dr.  Harris  Mayer),  1281  S.  We.stwood 

Blvd,  Los  Angeles  24,  Calif  Contract  AF  29(60I)-1897 

I  Republic  Aviation  Corporation,  MSD  (Mr.  Jeff  Meeker),  223  Jericho 

Turnpike,  Mineola,  Long  Island,  NY  Contract  AF  29(601)*2792 

I  Convair,  Division  of  General  Dynamics  (Dr.  Daniel  Hamlin).  3165 
Pacific  Highv/ay,  San  Diego  12,  Calif  Contract  AF  29(6011-2567 

1  University  of  Chicago,  Institute  of  Air  Weapons  Research,  Museum 

of  Science  and  Industry,  Chicago  37,  Ill,  THRU:  ASD  (W'.VRP), 

(Mr.  Card),  Wright- Patterson  AFB,  Ohio  Contract  A F33(i)  1  6) - 6824 

I  Armour  Research  Foundation  (Mr.  C.  Haaland),  10  W.  35th  St.. 

Chicago  16,  Ill 

1  Lincoln  Laboratories  (Dr.  J.  Panncll),  Massachusetts  Inst Hite  of 

Technology,  P.  O,  Box  73,  I.exington,  Mass 

1  Edgerton,  Germeshausen  f-  Grier,  Inc.,  (Dr.  M.  Shuler),  IfOO 

Brookline  Ave.  ,  Boston  15.  Mass 

1  Official  Record  Copy  (SWRPA,  l.,t  .Scliiff) 
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